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Abstract 

, We study the chiral two-matrix model with polynomial potential functions V and W, 

which was introduced by Akemann, Damgaard, Osborn and Splittorff. We show that the 
squared singular values of each of the individual matrices in this model form a determinan- 
tal point process with correlation kernel determined by a matrix-valued Riemann-Hilbert 
, problem. The size of the Riemann-Hilbert matrix depends on the degree of the potential 

function W (or V respectively). In this way we obtain the chiral analogue of a result of 
Kuijlaars-McLaughlin for the non-chiral two-matrix model. The Gaussian case corresponds 
' to V,W being linear. 

For the case where W{y) — ?/^/2 -I- ay is quadratic, we derive the large n-asymptotics 
<~i ' of the Riemann-Hilbert problem by means of the Deift-Zhou steepest descent method. This 

C^, proves universality in this case. An important ingredient in the analysis is a third-order 

(-H ' differential equation. 

Finally we show that if also V{x) — x is linear, then a multi-critical limit of the kernel 
exists which is described by a 4 x 4 matrix-valued Riemann-Hilbert problem associated to 
the Painleve H equation q"{x) — xq{x) + 2q'^{x) — v — 1/2. In this way we obtain the chiral 
analogue of a recent result by Duits and the second author. 

Keywords: chiral random matrix theory, two-matrix model, Riemann-Hilbert problem, 
' universality, multi-criticality, singular values, determinantal point process. 
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1 Introduction 

Chiral random matrix theory provides a useful calculational framework for solving various phys- 
ical problems, especially in quantum chromodynamics. Originally this was developed for the 
chiral 1-matrix model with a Gaussian potential [42 1 1441 H5] . 

Recently a chiral 2-matrix model was introduced by Akemann-Damgaard-Osborn-Splittorff 
[3]. The model is a Hermitian analogue of a non-Hermitian model introduced earlier by Osborn 
|40] . The motivation for this matrix model comes again from quantum chromodynamics. It 
turns out that the mixed correlation functions for n — )• oo strongly depend on a certain constant 
in the model, which has a physical interpretation as the 'pion decay constant' This suggests 
a way to compute F^^ via numerical lattice simulations. A numerical implementation of this idea 
is provided in |14j and some recent developments are, e.g., in [3B]. The method for calculating 
Fj^ was originally developed without random matrices \12\ 113] . The chiral two-matrix model 
yields an elegant theoretical framework that brings this method to its full power, see also [H [7] 
for the connection between the matrix model and the limiting physical theory. 

The chiral two-matrix model in [3] in its simplest, i.e. 'quenched' form, is defined by the 
probability distribution 

i exp (-nTr($*$ + ^'*^')) d$ d^, 

defined on pairs of rectangular complex matrices (<1>, ^'), both of size n x (n + ly), where n and 
v are integers, the superscript * stands for the conjugate transpose, and Z„ is a normalization 
constant. Here d<^ and d^' are the flat complex Lebesgue measures on the entries of ^ and ^. 
After the change of variables, see [4], 

^>1 = ^> + ^1^', ^)2 = ^> + ^2^, /il,/i2GK, 



2 



this becomes 



^ exp (-nTr(ci^'i$i + C2^<^2 - r($i$2 + ^2^1))) d$i d«>2, (1.1) 

for certain constants ci,C2 > and r € M and a normalization constant Z„. Here ci,C2,t are 
known functions of /xi and fi2- We assume without loss of generality that r > 0. For this 
measure to be finite we also need < < C1C2. 

As noted in [3], one can consider the more general model 

^ exp (-nTr(y($;^$i) + VF(«>2$2) - r($;^$2 + ^2^1))) d$i d$2, (1-2) 

where are polynomials with positive leading coefficient. This is the setting that we will 

consider in this paper. We assume without loss of generality that is a nonnegative integer. 

We note that the general model in [3] has extra determinantal factors in the probability 
distributions (|l.l|) ~ (|1.2p . What we consider here is the 'quenched' case where there are no 
determinants of this kind. In fact, the correlation functions of the general model can be expressed 
in terms of their quenched variants [H [6] . 

The chiral two-matrix model can be analyzed with the help of biorthogonal polynomials 
[11[26]. The monic polynomials are denoted Pj,n{x) and Qk,n{v)i of degree j and k respectively, 
and they satisfy the biorthogonality relations 

/ / Wn{x,y)Pj,n{x)Qk,n{y)'ixdy = Kk6j^k, Kk^O, j,k = 0,1,2,- ■■ , (1.3) 
Jo Jo 

with respect to the weight function 

Wn{x,y) = (:Ey)-/2l,(2Tn^)e-"(^(-)+^(j')), (1.4) 

see [H Eqs. (2.15) and (2.16)], with Ii, the modified Bessel function 

00 

The weight function is well-defined for all u > —1, not necessarily an integer. Following the 
approach by Ercolani-McLaughlin [25j one shows that the biorthogonal polynomials exist, are 
unique, and have real and simple zeros. In the case where V, W are both linear, the biorthogonal 
polynomials are explicitly known to be given by Laguerre polynomials [1]. 

Incidentally, there also exists a non-Hermitian version of the chiral two-matrix model, due 
to Osborn [401. In that case, Iu{x) in ()1.4p is replaced by the modified Bessel function of the 
second kind Ku(x), which is defined in (j4.12p below. 

Our interest lies in the singular values of the matrix $1. They form a determinantal point 
process with a correlation kernel called in [2]. We find it more convenient to study the 
determinantal point process of the squared singular values of $1 which is then described by the 
correlation kernel 

Kn{xi,X2) = ^ ( / Wn{xi,y)Qk,n{y)dyjPk,n{^2), (1-6) 

with Kk as in ()1.3p . There are three other kernel functions in [4j (as usual in the Eynard-Mehta 
setting [26]) but we will be only interested in Kn{xi, X2)- Indeed, in this paper we will express 
Kn{xi,X2) in terms of a Riemann-Hilbert (RH) problem. It is an open problem to apply our 
method to the mixed correlation functions, which describe the interaction between the singular 
values of $i and $2- 
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The kernel ()1.6p describes a well-defined determinantal point process on the positive real line 
for any value of > — 1. Our results on this point process will also hold for general v > —1 and 
have a random matrix interpretation in case v is integer, i.e. then the particles correspond to 
the squared singular values of the matrix $i in ()1.2p . In the following we will often refer to the 
process described by the kernel ()1.6p as the (squared) singular value process, even if u is not an 
integer. 

Summarizing, biorthogonal polynomials allow to obtain the correlation functions in the chiral 
two-matrix model. One can also find analytic expressions for the distributions of the individual 
singular values, see [21 [S]. 

In this paper we will show that the biorthogonal polynomials and Qk,n{y) can be 

characterized as multiple orthogonal polynomials \33\ H3] with respect to a suitable system of 
weight functions. Consequently we will express the kernel Kn{x,y) in terms of a RH problem. 
This yields the chiral analogue of a result of Kuijlaars-McLaughlin [33]; see also [8| 1251 [3T| for 
some RH problems of a different nature for the non-chiral two-matrix model. 

The paper [4J contains a detailed analysis for the case of linear potentials V{x) = cix, 
W{y) = C2y. In the present paper we focus on the quadratic case 



The other potential V{x) will be allowed to be an arbitrary polynomial with positive leading 
coefficient. Under these assumptions on V and W, we will be able to perform a Deift-Zhou 
asymptotic analysis of the RH problem, yielding the asymptotics of the kernel Kn- In this way 
we get the chiral analogues of the results by Duits-Kuijlaars-Mo |22l [231 [Ml [39] . 

Our results imply universality in the case of a quadratic potential W{y). Universality results 
for the chiral 1-matrix model were obtained in [21 IBD| [55] . 

Under the additional assumption that V{x) = x, we give an ar-phase diagram and discuss 
the phase transitions. In particular the phase diagram indicates a multicritical point for the 
values of parameters a = — 1 and r = 1. We describe the local behavior of the singular value 
process near this multicritical point by means of a triple scaling limit leading to the chiral version 
of the main result in [20] . The new critical kernel that we obtain in the scaling limit is expressed 
in terms of a 4 x 4 matrix- valued RH problem that was introduced by one of the authors in [16] 
to describe a critical phenomenon for non-intersecting squared Bessel paths. Our kernel will be 
built from the same RH problem, but in an essentially different way than in |16] . 

2 Statement of results 

Our first results hold in the general case where V and W are polynomial potentials with positive 
leading coefficients. We find it convenient to rewrite (jl.4p as 



W{y) = y^/2 + ay 



Wn{x,y) = fn{xy)e 



n{V(x)+W(y)) 



(2.1) 



with 



fn{x) : 



= X' 



,u 



'/2/^(2rnV^). 



(2.2) 



We also introduce the functions 




(2.3) 



for / G N U {0} and write 



r := degW — 1. 



(2.4) 
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2.1 Multiple orthogonality relations 

This is our first main result. 

Theorem 2.1. (Multiple orthogonality relations) 

(a) The hiorthogonal polynomial Pj^n{x) in ()1.3p is the type II multiple orthogonal polynomial 
with respect to the weight functions 



More precisely, 



e-"''(")Vn(x), / = 0,...,2r. 

P,>(x)x'=e-"^(^')V(x)dx = 0, 



(2.5) 



2r+l 



where we recall the notations (|2.3p - ()2.4p and [• 



for I = 0, . . . , 2r and k = 0, . . 
denotes the integer part of a number. 

(h) The polynomial Pj,n{x) also satisfies the alternative system of multiple orthogonality 
relations 

fOO 

Pj-n(x)x^'w;;,„(x)dx = 0, (2.6) 







for I = 0, . . . ,2r and k = 0, . . . 



2r+l 



with the weight functions 

e-"^(")/ii,„(x), / = 0,...,r, 

e-nV(x)^y^_^_^ ^(x), / = r + 1, . . . , 2r. 



(2.7) 



Theorem 12.11 will be proved in Section [3l A similar result holds of course for the biorthogonal 
polynomials Qj^n{x)- 

2.2 Riemann-Hilbert problem and correlation kernel 

Theorem 12.11 asserts that the polynomials Pj^n satisfy multiple orthogonality relations of type II. 
Hence, they are characterized by the following RH problem [43]. In the statement of the theorem 
we will use the system of multiple orthogonality relations ()2.6p rather than ()2.5p . We write 
M+ := [0,oo). 

RH problem 2.2. We look for a (2r + 2) x (2r + 2) matrix-valued function Y : C \ M+ — >• 
(j~<(2r+2)x(2r+2) gQ^igfying the following conditions. 

(1) Y{z) is analytic (entrywise) for z G C \ M^. 

(2) Y has limiting values Y± on M"*", where 1+ (Y^) denotes the limiting value from the upper 
(lower) half-plane, and these limiting values satisfy the jump relation 



y+(x) = y_(x) 



/I WQ^n{x) 

1 



W2T,n{x)\ 









(2. 



1 / 



for X G M+, where wi^n{x), I = 0, . . . ,2r, is given in (|2.7p . 
(3) As z — )• oo, we have that 

1 



Y{z) = [I + {-] ] diag(z",z-"«,z-"i,...,z-"2r)^ 



(2.9) 



where ni 



n+2r-l 
2r+l 



for Z = 0, . . . , 2r. 
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(4) Y{z) has the following behavior near the origin: 

Y{z)di8.g{l,h{z)-\...,h{z)-') = 0{l), 
Y-^{z)dmg{h-\z),l,. . . ,1) = 0{1), 

re the O 

denotes the inverse transpose, and 



as z — )■ 0, z G C \ M"^, where the O-condition is taken to be entrywise, the superscript 



zl", if-l<iy<0, 
h{z) = { log|z|, ifu = Q, (2.10) 
1, ifiy>0. 



This RH problem has a unique solution. It is constructed out of the type II multiple orthog- 
onal polynomials and their Cauchy transforms. In particular, the (1, 1) entry of Y is Pn,n- The 
inverse transpose Y~'^ is given in terms of the associated type I multiple orthogonal polynomials. 
We refer to [l3] for the details. 

Condition (4) is needed to ensure that the solution of the RH problem is unique. The exact 
form of (4) follows from an analysis of the formulas for Y and Y~'^ in terms of the multiple 
orthogonal polynomials mentioned above. This can be done as in [34|, Proof of Theorem 2.4]. 
An essential ingredient is that wi^n{z) = 0{z'^), as z — )■ for / = 0, . . . , 2r, which is immediate 
from ([22]) and (f2ni . 

The correlation kernel Kn in ()1.6p has the following representation in terms of Y, see 

1 T 

Kn{x,y) = ^^.^^ _ (O Wo^n{y) Wi^niy) ■■■ W2r,niy))^+^iy)Y+i^) {'^ ••• O) , 

(2.11) 

where the superscript denotes the transpose and both row vectors have length 2r + 2. 

This representation allows us to derive the large n limit of the correlation kernel in the case 
of a quadratic potential W{y) = y'^/2 + ay by performing a Deift/Zhou steepest descent analysis 
on RH problem 12. 2[ As will be clear from the analysis, this corresponds to a quartic potential in 
the non-chiral two-matrix model studied by Duits-Kuijlaars-Mo in \22\ 1231 [39]. We will largely 
follow the line of thought in these works, however, at several places complications will arise. 

The key to the steepest descent analysis is a third order differential equation for ho^n{x) from 
(|2.3p . to be stated in Section HI It plays the same role as the Pearcey equation in [221 [23] but is 
considerably more complicated. The steepest descent analysis itself will be described in detail 
in Section [5l 



2.3 Vector equilibrium problem and connection with the non-chiral two- 
matrix model 

From this point we will assume that the second potential is quadratic 

Wiy) = !j + ay, a G M, (2.12) 

and that V is a polynomial with a positive leading coefficient. As one of our main results we 
will characterize the limiting mean squared singular value distribution of $i in this case using 
a vector equilibrium problem. More precisely, it is the first measure of a triplet of measures 
minimizing an energy functional under certain conditions. Before we introduce this vector 
equilibrium problem we review a related equilibrium problem that arises in the study of the 
non-chiral two-matrix model. 
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Vector equilibrium problem for the non-chiral two-matrix model 

In the recent paper [23] . Duits, Kuijiaars, and Mo consider a Hermitian two-matrix model of 
the form 

exp (-nTr(y^C(Mi) + W^^{M2) - tMiA^)) dMi dMs, (2.13) 

n 

defined on pairs (Mi,M2) of n x n Hermitian matrices. Here, is a normalization constant, 
is a general even polynomial with a positive leading coefficient, W^'~' is a quartic polynomial 
given by 

w^''^(y) = ^/ + fy', aeM, (2.14) 

and r > is the coupling constant. Throughout this paper, we use the superscript ^'^ to 
distinguish functions and constants related to the non-chiral two-matrix model from similar 
notions in the chiral two-matrix model. 

We define, as usual (see [Hj), the logarithmic energy of a measure v by 

= /"/"log du{x) diy{y), (2.15) 

JJ \x-y\ 

and the mutual energy of two measures vi, U2 by 

I{yuy2) = //log di^2(y). (2.16) 
J J \x-y\ 

The main result of |23j is that the limiting mean eigenvalue distribution of Mi can be described 
by the first component of a triplet (/xj^^, yUg**^, /ig'*-'), which are three measures minimizing the 
energy functional 

3 2 

i?''^(^i,i^2,i^3) := ^/(^,) -Y,I{u,,v,^i) + j V^''{x)dvi{x) + j Vi'''{x)dv:,{x), (2.17) 

with V^^ and being certain symmetric external fields on M, and where the minimization 
is among all positive measures such that 

(a) ui is a measure on M with total mass 1; 

(b) U2 is a measure on iR with total mass 2/3 that satisfies the constraint 

V2 < a^^, (2.18) 
where (72'*^ is a certain positive symmetric measure on the imaginary axis; 

(c) 1/3 is a measure on M with total mass 1/3; 

(d) < 00 for j = 1,2,3. 

This equilibrium problem clearly depends on the input data V^'^, V-^^, and cr2^^. For the exact 
definitions of these notions we refer to [23]. The unique solvability of the vector equilibrium 
problem follows from [23\ [29] . 
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Vector equilibrium problem for the chiral two-matrix model 

Given the triplet {V, W, r) defining the chiral two-matrix model, we define 

V^C{x) = ^V{x^) and W^^ (y) = (y^) . (2.19) 

Note that this definition is consistent with ([231]) and (f2J2]) . The triplet (y^^,VF^^,r) then 
characterizes an associated non-chiral two-matrix model. 

The vector equilibrium problem that is appropriate in our chiral setting is a 'squared' version 
of the above vector equilibrium problem for the associated non-chiral two-matrix model. More 
precisely, we consider the energy functional 

E{ui,iy2,U3) ■■=^H''j) - ^H^j^'^j+i) + J Viix)dui{x) + J V3{x)du3{x), (2.20) 
j=l j=i 

and define the input data 

Vi{x):=2V,''''{V^), Vs{x):=2Vi'''{V^), da2(,x) := 2 da^"" (iV^) . (2.21) 

Then Vi and V3 are external fields on M"*" and a2 is a positive measure on M". The vector 
equilibrium problem is then to minimize E[ui, 1^2, 1^3) among all positive measures z^i, z/2 and z/3 
satisfying the following conditions. 

(a) z/i is a measure on with total mass 1. 

(b) z/2 is a measure on := (—00, 0] with total mass 2/3 such that 

z/2 < CT2, (2.22) 

where o"2 is defined in ()2.2ip . 

(c) z/3 is a measure on with total mass 1/3. 

(d) /(z/j) < 00 for j = 1,2,3. 

This vector equilibrium problem has a unique solution described in the following theorem. 
We denote the support of a measure z/ by 5(z/). 

Theorem 2.3. The equilibrium problem ()2.20p - ()2.22p has a unique solution {fii, 112, f^s). More- 
over, if ^ S{fii) or ^ S{a2 — ^2) then 

N 

S{l^i) = \J[aj,bj], 
i=i 

for some N £ N and < ai < 61 < 02 < • • • < Qat < bj\f and on each of the intervals [aj, bj] in 
S{ni) there is a density 

dpi _ , . _ \l9jix)^ibj - x){x- ttj), x G [aj,6j], if Oj > 0, 

dx "^'^^^ " \ lg^(x)y/{bi-x)x-\ XG[0,6i], if O^ = 0, 

where gj is nonnegative and real analytic on [aj,bj]. 
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Proof. We claim that 

(d^i(x), d^2(x), d^3(x)) := (2d/ifC(V^),2d;U^C(V^),2d;u|'C(V^)), (2.23) 

where the right hand side denotes the solution to the equilibrium problem ()2.17|) - ()2.18p with the 
same parameters r, a. To prove (j2.23p . let us denote for any measure /x symmetric with respect 
to the origin, the squared measure /2 by the rule d/2(x) = 2dfj,{y/x). Since fi is symmetric we do 
not have to precise the branch cut of the square root. Thus if fi has a density dfi{x) = p{x) 6.x 
then 

d/I(x) := 2p{^) dV^ = dx. 

yx 

Now the mutual energies of two measures vi , U2 and the corresponding squared measures vi , V2 
are related by 

I {1^1,^2) = 2/(1^1, 1/2), 

see e.g. [HI Th. IV.l. 10(f)]. With the help of this relation, and using fn\\i . (|2T4]) and (|2T^ . 
the claimed relation ()2.23p follows. 

Given this, the theorem is a corollary of [23^ Theorem 1.1]. □ 

Remark 2.4. If ai = the density of /ii blows up like an inverse square root at the origin as is 
the case for the Marchenko-Pastur density [37] . 



2.4 Classification into cases 

In this model we distinguish a number of regular and singular cases depending on the supports 
and densities of the measures /xi, (72-/^2, and /is. The classification in our chiral setting 
is inherited from the classification for the associated non-chiral models given in \12>\ Section 
1.5], e.g. we say that our chiral model belongs to Case I if the associated non-chiral model 
is in Case I according to [23l Section 1.5]. This leads us to distinguish five generic cases and 
8 singular cases. The classification of the five generic phases depends on whether is in the 
support of the measures /^i, o"2 — /X2i and /X3 or not. We have the following cases with regular 
behavior of the supports at zero: 

Case I G 5(/xi), ^ S(a2 - 1^2) and G ^(/is), 

Case II ^ S{ni), ^ S{a2 - ^2) and G S{n3), 

Case III ^ S{fii), G S{a2 - H2) and ^ ^(/ig), 

Case IV G 5(/xi), ^ S{a2 - /X2) and ^ 8(113), 

Case V ^ 5(/xi), ^ S{a2 - [12] and ^ ^(/ig). 

We discuss only three of the singular cases. For the remaining five singular cases we refer 
to [23]. The critical phenomena corresponding to these remaining five cases can already be found 
in the one-matrix model, so we will not discuss them here. 

Singular supports I G S{fii) fl S{a2 — /X2), ^ 3(^3). 

Singular supports II ^ 5'(/ii), G S{a2 - p.2) n S{ii3). 

Singular supports III G S{pi) n S{a2 - P2) n ^(/is). 

The last one is a multi-singular case. 

Except for Case V, all these phenomena already occur in the simplest version of the model in 
which the potential V{x) = x. For this situation we will establish a phase diagram in Section[2]6j 

The next theorem will only be proved in the regular cases that we define as follows. 
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Definition 2.5. The triplet {V, W, r) is regular if the associated triplet (y , W ,t) is regular 
in the sense of [231 Definition 1.3], with V^'^ and W^'^ as in (l2J9l) . 

In particular for a regular triplet it holds that S[fii)r\S[a2 — ^12) = and S'(//3)nS'((T2 — ^2) = 
0. Moreover the functions gj, j = 1,. . . ,N, in Theorem 12.31 are nonzero for all x S [aj,bj], i.e. 
the density of the measure /xi does not vanish in the interior of the support and behaves like a 
square root at the nonzero endpoints of the support. The same holds for the densities of a2 — fJ'2 
and /i3. There is an extra condition on the variational inequality for which guarantees that 
no extra interval emerges in the support of ^ui when continuously varying the potentials. 



2.5 Limiting mean singular value distribution 

The measure /ii is the limiting mean squared singular value distribution of the matrix <l>i in the 
chiral two- matrix model as n — )• 00. This statement holds for W as in ()2.12p and general V and 
T, but we will only prove it for the regular cases, see Definition 12.51 

Theorem 2.6. Suppose {V, W, r) is regular. Let fii be the first component of the minimizer 
(/xi, /X2, /^s) of the vector equilibrium problem (|2.2U|) - (|2.22|) . Then /ii is the limiting mean particle 
distribution of the determinantal point process with correlation kernel ()1.6p as n —)• 00 with n = 
mod 3. 

This theorem is equivalent to the statement 

lim —Kn{x,x) = pi{x) = — ^(x), X > 0, n = mod 3. 
n^oo n ax 

The proof is given in Section [5^ It is based on a lengthy Deift/Zhou steepest descent analysis 
carried out in Section^ Without too much extra effort one could also obtain the universal scaling 
limits that are typical for unitary random matrix ensembles. More precisely, if pi{x*) > 0, where 
pi denotes the density of the measure pi, we retrieve the sine kernel as a scaling limit 

1. . * - * i7 . sin7r(x — y) 

limi^„(x* + r—:,x*+ - i- \ yj 



npi{x*)^ npi{x*) J Tr{x — y) 

If X* is a nonzero endpoint of S{pi), i.e. x* G {oi, 61, . . . ,07V, b^} \ {0}, then 

r 1 f *^ ^ *, y \ Ai(x)Ai'(y)-Ai'(x)Ai(y) 

lim - — —j-Kn I X ± - — 7^77^)2; ± — 



where we choose the +-sign ( — sign) if x* is a left (right) endpoint of the support of pi. Recall 
that we are in the regular case so that the density vanishes like a square root at x* > 0. However, 
if x* = ai = then the density blows up as an inverse square root and we would obtain the 
Bessel kernel of order as a scaling limit 

lim -^—K ( ^ y \ - - V^J'v{\^)Jy{y/y) 



n^co {^cuY \{cn)'^^ (cn)"^ J 2{x — y) 

for x,y > Q and a suitable constant c > 0. We will not discuss this any further. 



2.6 Phase diagram in the quadratic/linear case 

Prom here we restrict ourselves to the very specific model of quadratic and linear potentials 

V{x) = x, W{y) = ^ + ay, a e M. (2.24) 

For this concrete model we can construct a phase diagram, i.e. determine which values of (q, t) 
correspond to which case of the classification in Section ^IM It turns out that the case 'Singular 
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supports III' occurs and this will get most of our attention. To establish the phase diagram, 
we first discuss the behavior of the supports of the measures /Ui, ct2 — jU2, and fis and how they 
depend on a and r. It follows from ()2.23p and that the supports of the measures fii, o"2 — /X2, 
and //3 have the following form 

supp(/ii) = [/3i,/3o], 

SUpp((T2 - /U2) = (-00, -^2], 
SUpp(/i3) = [/33,00), 

for some /3o > /3i > 0, /32,/33 > that all depend on the values of q € M and r > 0. We 
distinguish a number of cases, depending on whether /3i, /32, or are equal to zero, or not. At 
least one of these is zero, and generically, no two consecutive ones are zero. According to the 
classification in Section 12.41 we have 

Case I: /3i = 0, /32 > 0, and jS^ = 0. 

Case II: pi > 0, /32 > 0, and ^3 = 0. 

Case III: pi > 0, /32 = 0, and /S^ > 0. 

Case IV: pi = 0, /^a > 0, and P3 > 0. 

Case V does not occur in this specific model. 

The phase diagram in Figure [1] shows which values of (a,T) correspond to these four cases. 
The different cases are separated by the curves given by the equations 



r = \/a + 2, —2<a< 00, and r = \ , —00 < a < 0. 

V a 

On these critical curves two of the numbers Pi,P2 and P^ are equal to zero. For example, on 
the curve between Case III and Case IV, we have Pi = P2 = 0, while P3 > 0. Finally, note the 
multi-critical point (a, t) = (—1,1) in the phase diagram, where Pi = P2 = Ps = 0. All four 
cases come together at this point in the (a, r)-plane. The nature of this multi-critical point is 
discussed in the next section. 

We do not study the other types of critical behavior in this paper, but we shortly comment 
on them now. The geometry of the supports of fii, cr2 — Ai2i and fi^ suggests that the transition 
on the curve r = \/a + 2, a > —1, is described by the inhomogeneous Painleve II kernel as in 
|10j . The transition on the curve r = y/—l/a, a < —1, on the other hand is described by the 
chiral version of the Pearcey kernel as in [52] . The transitions on the dashed lines in Figure [T] 
are situated on the non-physical sheets of the underlying Riemann surface and, hence, do not 
affect the local correlations of the singular values of 

The curve r = \/a + 2, a > —2 is also critical for the singular values of ^2- The other curve 
is not critical in that context. 

Figure [2] shows the phase diagram again together with the shape of the limiting mean squared 
singular value densities above, below, and on the critical curve. 

2.7 A triple scaling limit 

Let us now focus on the squared singular value process of $1 near the critical parameters r = 1 
and a = —1, by means of a triple scaling limit. To this end, we rescale a and r near the critical 
values in the following way 
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Figure 1: The phase diagram in the ar-plane: the critical curves r = \/a + 2 and r = y — ^ 
separate the four cases. The cases are distinguished by the fact whether is in the support of 
the measures fJ-i, 0-2 — H2, and ^3, or not. 

for a, 6 G M. We also scale the space variables with 

x = un~^^^, and y = vn~'^^^, u,v > 

and compute the limiting behavior of Kn{x,y) as n — )• 00. 

The limiting kernel is characterized by the solution to the following RH problem introduced 
in [16]. The RH problem has jumps on a contour in the complex plane consisting of 10 rays 
emanating from the origin. More precisely, we fix two numbers '^i , (p2 such that < (pi < (p2 < 
tt/2 and define the half-lines F^, A; = 0, . . . , 9, by 

Fo = M+, Fi = e*'^iM+, F2 = e''^m+, F3 = e''^^-'P^^R+ , F4 = e'^''-'^^'>R+ , (2.26) 

and 

F5+fc = -Ffc, fc = 0,...,4. (2.27) 

All rays F^, A; = 0, . . . , 9, are oriented towards infinity, as shown in Figure [3l We also denote by 
r^fc the region in C which lies between the rays F^ and F^+i, for k = 0, . . . ,9, where we identify 
Fio := Fq. Now we can state the RH problem. 

RH problem 2.7. We look for a 4 x 4 matrix-valued function M : C \ (^Ufe=o) ~^ C^^^ 
(which also depends parametrically on v > —1/2 and on the complex parameters ri,r2,s,t £ C) 
satisfying the following conditions. 



(1) M(C) is analytic (entrywise) for G C \ ^Ufc=o-'^fc) • 

(2) For C G Ffc, the limiting values 



M+(C)= lim M(z), M_(C)= lim M{z) 

z on +-side ofF^. z on —-side o/F^ 

exist, where the +-side and — -side ofT^ are the sides which lie on the left and right ofT^, 
respectively, when traversing Tk according to its orientation. These limiting values satisfy 
the jump relation 

M+(C) = M_(C)Jfc(C), A; = 0,...,9, 
where the jump matrix Jk{C) for each ray F^ is shown in Figure\^ 
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r 




Figure 2: The phase diagram in the ar-plane with density plots. Below the critical curve the 
limiting mean squared singular value density of is supported on one interval touching the 
origin. At the origin the density blows up like an inverse square root. Above this line the singular 
values cluster on an interval away from zero. On the critical curve the transition between these 
two regimes occurs. There the density is supported on an interval including zero. On the part 
of critical curve to the right of the point (—1, 1) (dotted) the density vanishes like a square root. 
To the left of this point (dashed) the exponent at the origin is —1/3. At the point (—1, 1) we 
have a transition between these two types of critical behavior. There the density blows up with 
an exponent —1/4. 



(3) As C, ^ oo we have 

M(C): 



Ml M2 



diag((-c)-^/^^^/^(-c)^/^c^/^) 



X ^diag (e-^2(C)+*C,e-V'i(C)-*C^eV2(C)+*C^eVi(C)-*C^ 



where the coefficient matrices Mi,M2, . . . depend on the parameters v, ri,r2, s and t, but 
not on C, and where we define 



V2 



( \ -i 0\ 

1 i 

-i 1 

Vo z ly 



(2.28) 



V'i(C) = ^riC^/^ + 2<V2^ = ^r2(-C)'/' + 2s(-C)^/'. (2.29) 
(J^) M{Q behaves /or C ^ as 

M(C) = 0(n, M-\0 = OiC), ^/P<0, 

and 

' M(C) diag(C-^ C"') = 0(1), C G ^^1 U Og, 

< M(C)diag(C,r",r",n = C?(l), CeOgUOe, i/i>>0. 

. M(C) = 0{Cr% C ^ (J^i U U ile U J^s), 



13 







Vo 



) 


0\ 












) 1 







1 


l) 






A 









1 







1 




^0 


e'^'^* 



Figure 3: The figure shows the jump contours in the complex i^-plane and the corresponding 
jump matrix on F^, A: = 0, . . . , 9, in the RH problem for M = M{(). We denote by 0^ the 
region between the rays F^ and T^^i. 



It was proved in [TB] that this RH problem is solvable for ri = r2 > and s,t S M. In 
the same paper this RH problem is shown to be related to the Hastings-McLeod solution of the 
inhomogeneous Painleve II equation q"{x) = xq{x) + 2q^{x) — u — 1/2. 

We transform the RH matrix M{z) into a new matrix M{z) as follows 

M{z):=d\B.g{z^l\z-y\z^l\z-^l')dmgi^(^^ ' (-1 l)) ^(^^'^')- (2.30) 

The transformed matrix M{z) depends on the same parameters ri,r2, 3,1,1/ as M(z). The 
matrix M{z) satisfies a RH problem by itself but we will not state it here. 
For u,v > 0, we now define the critical kernel K^'^{u, v; s, t, u) by 

K^''(u,v;s,t,u) = -(0 -ie-"^' l) M (v; s,t,i)r^ 

2m[u — V) 

xM{u;s,t,i)){0 -ie"''' l)^, (2.31) 

where M{z; s, t, 9) is defined by ()2.30p with M{C,) = M(C; s, t, u) the unique solution to RH 
problem 12.71 with parameters 

n = r2 = 2, 

< s,t E M, (2.32) 
u = v + l/2. 

It can be shown as in [20l Appendix A] that this kernel is different from the hard-edge tacnode 
kernel discovered in |16j . 

We can then state our final main result. 

Theorem 2.8. Assume v > —1 and let Kn be the kernel in (jl.6|) with V and W as in (j2.24p . 

If I' is integer this kernel describes the squared singular values of Mi when averaged over M2. 
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Set 

for a, 6 G M. Then for n — >• oo and n = mod 3 

uniformly for u,v in compact subsets ofM^. 

This theorem will be proved in Section 16.81 The prefactor v^l'^ju^l'^ has no influence on 
the singular value correlations. The phase diagram in Figure [1] suggests that by taking proper 
scaling limits, this critical kernel converges to the inhomogeneous Painleve II kernel or to the 
kernel in [32]. The non-chiral version of this statement has been shown in [20^ I28j . 

2.8 Organization of the paper 

The rest of this paper is organized as follows. In Section [3] we prove Theorem l2.1l on the multiple 
orthogonality relations for the biorthogonal polynomials and Qk,niy)- In Section H] we 

obtain and analyze a third-order differential equation related to the chiral two-matrix model in 
the case of quadratic W{y) = y^/2 -|- ay. In Section [5] we perform a steepest descent analysis on 
RH problem 1 2. 2 1 for quadratic W leading to the proof of Theorem 12. 6i In Section [6] we adapt the 
steepest descent analysis of the preceding section to prove Theorem 12.81 The main difference is 
in the local parametrix at the origin which is now built using the solution to RH problem 12. 7[ 

3 Proof of Theorem [231 

3.1 Preliminary lemma 

Lemma 3.1. The function fn{x) in (|2.2|) satisfies the differential equation 

xf:{x) -{u- = (rnffnix). (3.1) 

Proof. This is immediate from the definition ()2.2|) and the fact that luiz) satisfies the modified 
Bessel equation [1] 

z^I'liz) + ziliz) - {z^ + u^)I,{z) = 0. (3.2) 

□ 

3.2 Proof of Theorem [2JTa) 

We claim that for any polynomial P and any nonnegative integer /, one has 

/ x'=P(x)e-"^(^) / y'fn{xy)e-^'^^y^dydx = / / P(x)7rfc,K2/)^n(x, y) dx dy, (3.3) 
Jo Jo Jo Jo 

where iTk^i is a polynomial of exact degree k{2r + 1) + /. 

We prove the claim by induction on k. The case A; = is trivial. Let us assume that the 
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statement is valid up to k and compute the left-hand side of (|3.3p for A: + 1 



1 f'OO f'OO 



u-l 



(rn)2 

-2 / x^P(x)e-"^(^) / 7fi+i+2r(y)/n(xy)e-"^(^)dydx 
ny Jo Jo 



where tTj denotes a polynomial of degree We used ()3.ip for the first equality. The second 
equality is based on integration by parts. Note that the integrated terms vanish for > and 
cancel each other out for —1 < z> < 0. Now expanding 



1 u-l l+l+2r 



i=0 

and applying the induction hypothesis yields 



OO poo 

x^+ip(x)e-"^(^) / y^fn{xy)e-''^^yHydx 







l+l+2r 



' roo poo 

T^aJ x^P{x)e-^''^-n y7n(xy)e-"'^(^)dydx 

-l+ 2r r-oo 

yZ / P{x)TTk,i{y)wn{x,y)dydx 

OO poo 

I Pix)7Tk+i,i{y)wn{x,y)dydx. 
Jo 

Here TTk+i,i is defined as 

l+l+2r 
1=0 

and is, therefore, of degree (A: + l)(2r + 1) + / as it should be. This proves the claim. 

To show ()2.5p . note that the set A = {TTk^i | A; = 0, 1, 2, 3, ...,/ = 0, ... , 2r} spans the space 
of polynomials. If we take P{x) = Pj^n{x) then the right-hand side of (|3.3p vanishes for any 
TTfc^/ G A satisfying k{2r + 1) + I < j. This leads to vanishing of the left-hand side of (|3.3p for 

every A; = 0, . . . , ^2r+i ' ^ ~ 0, . . . 2r. Hence, Pj^n is the j-th multiple orthogonal polynomial 
satisfying (|2.5p . This proves Theorem 12.1 I f a). 

3.3 Proof of Theorem ICTb) 

To prove ()2.6p . it is sufficient to establish linear relations between the functions hj^n and their 
derivatives. 

Lemma 3.2. The function xh'i^{x), see ()2.3p for hi^n, can he written as a linear combination 
of the functions hi^jn{x), j = 0, . . . ,r + 1, 

r+l 

x/i; = -(/ + l)/ii,„(x) + ^^ ^ Cjhi+j^n{x), (3.4) 
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for any / G NU {0}, where we write the polynomial W in ()2.ip as 

r+l 

W{y) = J2-y', (3-5) 

with coefficients Ck G K, c^+i > 0. Here we assume without loss of generality that the constant 
term of W vanishes. 

Proof. A straightforward calculation gives 

xhl^ix) :=xj y'_(/„(xy))e-"'^(2') dy 

poo Q 

= / y'^'^ifnixy))e--^(y^dy 
Jo dy 

= -^°°/n(xy)|(y'+V-^(^)) dy 

r+l 

= -(/ + l)/lz,„(x) + n^Cj/l;+j-„. 

i=i 

Here the first step follows from the definition ()2.3p . the second step is obvious by symmetry, the 
third step is integration by parts, and the last step uses (j3.5p . □ 

The above lemma then implies the alternative system of weight functions for the multiple 
orthogonal polynomials Pj^nix) in ()2.6p . This proves Theorem 12. If b). 

4 Differential equation for ho^n{x) with quadratic W 

From here we assume again that W is quadratic and depends on the real parameter a as in 
(I2H. Then r = degVF - 1 = 1. 

In this section we obtain a third-order differential equation for the function 

ho.n{x)= / /„(xy)e-"^(^)dy, (4.1) 
Jo 

see (j2.3p and (j2.2p . We will study this differential equation in detail and construct a Wronskian 
matrix from three special solutions to the equation. We will show that this Wronskian matrix 
satisfies a certain 3x3 RH problem, which should be regarded as a chiral analogue of the 3x3 
RH problem associated with the Pearcey differential equation used in [22^ [23] . 

4.1 Third-order differential equation 

In this section we prove the following result. 

Proposition 4.1. (Third order differential equation) Let W be quadratic as in ()2.12p . Then 
the function p{x) := hQ n{x) in (j4.ip satisfies the third order differential equation 

x^p" (yx) + (2 — 2u)xp"{x) + {an^T^x + z/^ — u)p'{x) — {r'^n^x + T'^n^va)p{x) = 0. (4.2) 

We note that a similar differential equation appears in a Mehler-Heine type formula in [19]. 
To prove Proposition 14. H we start with two lemmas. 

Lemma 4.2. The function ()2.3p can he expressed in terms of hi^n 

{mfh^i^r^ix) = xh'l^ix) -{u- l)h[^^{x), (4.3) 

for any I G NU {0}. 
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Proof. This is obvious from ()3.ip and ()2.3p . □ 

Lemma 4.3. Let W be quadratic as in (|2.12p . Then the functions hi^n, ^ = 0, 1, 2, 3, satisfy the 
relation 

/i3,„(x) + 2ah2,n{x) + {a^ - {v + 2)n-^)hi^n{x) - {{v + l)n-^a + t2x)/io,„(x) = 0. (4.4) 
Proof. We first assume that v > Q and calculate 

/oo 
y ((y + af - /„(xy)e-"^(s^) dy 

^2 







y yfn{xy)f^ (^e~-W(y)-^ (4.5) 

/>oo 

/ X (xy/;'(xy) + 2/;(xy)) e'^^^^) dy, 



where the first step follows from the definition (|2.3p , the second step uses (j2.12p , and the third 
step is integration by parts. Similarly, 

/ii,„(x) + a/io,n(x) = / {y + a)/„(xy)e-"^(^) dy 





-n-' fn{xy)^^ (e-"^(^)) dy (4.6) 



n 



dy 

W xf\xy)e-^'^^yHy. 



oo 







The lemma then follows from ()3.ip by a little calculation. Finally, note that the assumption 
V > Q was needed to make sure that the integration by parts in (j4.5p - ()4.6p creates no integrated 
terms. However, if < 0, integrated terms are present, but cancel each other out so that ()4.4p 
remains valid. □ 

Now we are ready to prove Proposition 14. 1[ 

Proof of Proposition First observe that ()3.4p and ()2.12p yield the relations 

x^Q^nix) = -ho^nix) + na/ii,„(x) + n/i2,n(x), (4.7) 
xh'i^nix) = -2hi^n{x) + nah2,nix) + n/i3,„(x), (4.8) 

and so on. Multiplying (|4.7p with an"^ and ()4.8p with n~^, and adding them up, we get 

h3,n{x) + 2a/i2,n(x) + Q^/ii,„(x) = (ax/i'o „(x) + Q/io,n(x) + xh'i „^{x) + 2hi^n{x)) . 

Substituting this into (|4.4p , we find 

xh'^ „(x) — vhi^n{x) + ax/i'o ^(x) — [va + T'^nx)ho^n{x) = 0. (4.9) 

On account of (|4.3p with / = this yields 

(Tn)^/ii,„(x) = x/io„(x) -{u- l)/io„(x), 
(Tn)2/i; „(x) = x/io'„(x) - (z^ - 2)/i'o',„(x). 

Inserting this in ()4.9p we get 

x(Tn)"2(2./jW^(2.) - (i/ - 2)/io„(x)) - z^(rn)-2(x/io^„(x) - (z^ - l)/i'o„(x)) 

+ ax/i'o „(x) — (z^a + T^nx)/io,n(2;) = 0, 
which with the notation p{x) := /io,n(a;) is equivalent to (|4.2p . This proves Proposition 14.11 □ 
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Remark 4.4. A little calculation shows that for any solution p{x) to ()4.2p the transformed 
function 

q{x) := x^^-^^'^/^p{x^) 

satisfies the differential equation 



3 //// \ , ( A 2 2 3 (^'^ + ^)'^ \ 1 1 \ 
X q [X) + ( 4r n ax x I q [x) 



16Z^3 + 60zy2 ^ 121/ - 7 ^ 4 3 4 4 2 2 2,n , ^^\ r ^ n 

For the special case v = —1/2 this reduces to the third order differential equation 

q"{x) + Ar'^n'^aq' (x) = 8T^n^xq{x), 

which appeared before in [231 eqn. (5.3)]. In that paper this ODE is solved using Pearcey 
integrals. 

4.2 Three special solutions to (14. 2 p 

We now construct three contour integral solutions po, pi, and p2 to the third-order differential 
equation ()4.2p . In the next section we will study the asymptotics of these functions as z — >• oo. 
The function po is defined as 

Po{z) := ho,n{^) = / /„(zy)e-"^(^) dy, (4.10) 







where we recall the definition of fn (j2.2p . Then pQ is analytic in C \ . 
The second solution pi to (|4.2p is defined as 

• i-oo 

p,{z):=- (zy)'^/2^,(2rn^)e-"^(^)dy, (4.11) 

J-oo 

where 

K,^,) = y-A^lfl^, zGC\(-oo,0] (4.12) 
2 sin(z^7rj 

is the modified Bessel function of the second kind. The right hand side of (|4.12p is replaced by 
its limiting value if u is an integer or zero. Then pi is analytic in C \ M. 
Finally, the third solution to (|4.2p is defined as 



P2{z) := r fn{zy)e-^'^^y^ dy. (4.13) 

J — oo 

This function is analytic in C \ M^. 

The functions po,Pi,P2 ah satisfy the differential equation (|4.2p . For po this follows di- 
rectly from Proposition 14.11 The same proof works for p2- Finally, for pi we first observe that 
{xY^'^ K^{2Tn^/x) satisfies the same differential equation ()3.ip as fni^) since also Ki, solves the 
modified Bessel equation (j3.2p . This then allows the proof of Proposition 14. 1 1 to be applied word 
for word. 

Lemma 4.5. The functions Po,pi,P2 satisfy the jump relations 

Po,+ (x) = e2'^>o,-(a;), XGM-, 

pi,+ (x) =pi-{x)-e-'''''p2-{x), X £ M+, 

Pi,+ {x) = pi,-{x) + e''>o,-(a;), X e R~, 

P2,+ ix) = e-2'^>2,-(2;), X E M+, 

where the subscripts + and — denote the boundary values obtained from the upper or lower half 
plane of C, respectively. 
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Proof. This follows from the definitions ()4.10p , ()4.1ip , and ()4.13p ; the connection formula ()4.12p ; 
and the fact that Iu+{x) = /j^_(x)e^'^'^* for x < which is immediate from (jl.Sp . □ 

4.3 Asymptotics for the solutions to f l4.2p 

Our next goal is to find the large z asymptotics of po, pi, and p2- In analogy with [23j we first 
introduce three functions 6j{z), j = 1,2,3, that will appear in these asymptotics. 
First we define the constants x*{a) and y*{a) that depend on a and r. We put 

x*{a) = { °' 4 ...3 " J n = -x%-a). (4.14) 

So if Q > we have that x*{a) = and y*{a) < whereas in the case a < it holds that 
x*(q) > and y*{a) = 0. 

Lemma 4.6. For every a € M and r > there exist analytic functions 6j : C\M — )• C, j = 1, 2, 3, 
with the following properties. 

(a) The jumps of the functions 9j are taken together in terms of the jumps for the diagonal 
matrix 

e{z) = diag{e,{z),e2{z),e3{z)), zec\M. (4.15) 

We have 

e+(x) 

< G+(x) 
G+(x) 

with x*{ol) and y*{a) as in ()4.14p . 

(b) We have as z — t- oo within 

2 3 

/or a constant D G M and j = 1,2,3. i7ere lo := exp(27ri/3). T/ie behavior in follows 
from the relation 9j{z) = Oj{z), j = 1,2,3. 

Proof. Define 

%(z) = 2^fC(^)^ j = 1,2,3, (4.16) 

where the 6^'~^ refer to the ^-functions used in [23j with potentials given in ()2.19p . From (2.5) 
and (2.8) in [23], it is readily seen that 

This, together with \1?>\ Corollary 2.2 and Lemma 2.4], implies our lemma. □ 

We are now ready to investigate the large z asymptotics of pq, pi, and p2. Let us first give 
a heuristic argument. 
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Lemma 4.7. For any solution p{z) to ()4.2p . there exist j G {0, 1,2} and C E C \ {0} such that 



p{z) = Cz^e"''^+^(^) ( 1 - -^u:^z-y^ + Du^h-^l'' + O(z-i)) , (4.18) 

as z ^ oo in closed sectors of the first quadrant {2; G C | Re > 0, Imz > 0}. Here D £ M. is a 
constant, and Oj is defined in Lemma \4-t>\ 

Proof. Substituting p{z) = e"^(^) in (j4.2p we get the following nonlinear ODE for f{z) := F'{z) 

{z^f\z) - T^z + ar'^zfiz)) + n'^ {2,z^ f [z) f {z) - 2{u - l)zf\z) - vr^a) 

+ _ 2{v - l)zf'iz) + uiu - l)f{z)) = 0. 

This ODE has solutions / with expansions 

-2^/ / 2\ 1 -1/3 2i -2/3 , - 1 _i a(9i/ + na^) .^/o 
' ^ 3 3n 81n 

+ {27u{u + 1) - 9 + 9(z. + l)a^n + qV) z'^/^ ^ Qiz'^), 

as z — )• 00, for any j S {0, 1, 2}. If Imz > 0, after integration this yields 

F{Z) = e,+^{z) + ^ + ^ log C- - + I)'^2,^-2/3 + 0(^-1)^ 

where C is the integration constant and D' € M. Then = e"^^^-* satisfies ()4.18p . □ 

Next we specialize Lemma 14.71 to the functions poi Pi^ and p2- 
Lemma 4.8. For each j = 0, 1, 2, we have 

pj,) = _^^^e"^-/3^2,^-^gn^^,+,(.) f 1 _ ^^^,^-1/3 + ^^2,^-2/3 ^ ^(^-1) 

V3n ^ 3r2/3n 

as z —)• oo in closed sectors of the upper half plane {z G C | Imz > 0}, where 

Co = l, C7i = -ze"f(2-+i), C2 = e-f(2-+i), (4.19) 

and is the constant in ()4.18p . /n i/ie lower half plane we have the same expansions hut with 
OJ replaced by u:~^ and Co,Ci,C2 by Co,—Ci,C2 respectively. 

Proof. First assume t = n = 1. Consider the ODE 

zq"'{z) + vq"{z) + aq{z) - q{z) = 0. (4.20) 
As observed in [19j, for any solution q{z) of ()4.20p 

p{z) = z^q'iz) 

solves (|4.2p . In [52] solutions to (|4.2Up are studied. The authors consider four solutions of the 
form 

qj{z)= f rt^-T+-*dt, i = 1,2,3,4, (4.21) 

where the contours Tj are shown in Figure HI 

The branch cuts for t'^~'^ are specified as follows: for j = 1 we take — 7r/2 < argt < 7r/2, 
for j = 2 we have 7r/2 < argt < 37r/2, for j = 3 we choose < argt < vr, and for j = 4 we 
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r4 



Figure 4: The contours of integration Tj, j = 1,2,3,4, used in the definitions (j4.2ip of the 
functions qj. 



put — vr < argt < 0. Note that our definition differs from the one in [32] by a multiphcative 
constant. 

We have the following relations 

q2{z)=q3iz)-e^'''''q4iz), for Rez > 0, (4.22) 

q^{z) = q3{z) - qiiz), for Rez < 0. (4.23) 

It is then clear that the functions 

z-'q'^iz) = z'' t-ie^-^+^* dt, J = 1, 2, 3, 4, 

generate the three dimensional solution space of (j4.2p . Moreover [32] provides detailed asymp- 
totics of these functions, e.g. as z — )■ oo such that < aigz < 7r/4 we have 

z'^q'liz) = ^/2^e-"V6_^^2^,f .^/^-..V3 ^ ^(^-i/S)^ ^ (4,24) 
z'^qUz) = ^/2^e-"V6e^_^^,^ef-^/^— ^^^/^ + 0(^-i/3)) , (4.25) 
z'^q'iiz) = -V2^e-"^/6e-^^a;2^^^|.^.^/3-a..V3 ^ o{z-'/')) . (4.26) 

The asymptotic behavior of z'^q2{z) follows using the connection formula ()4.22p . The idea 
is now to express po, pi, and p2 as linear combinations of the functions z'^q'-[z) to obtain the 
asymptotic behavior. This is done in Lemma 14.91 below. It leads to the following results 

p,{z) = -i.e"V3^^^f .^/3-..V3 ^ o{z''/'^)) , (4.27) 

as 2 —7- 00 such that < argz < 7r/4. This proves the lemma for r = n = 1 and in the sector 
< argz < 7r/4. The results in the remaining sectors are proved similarly. The lemma for 
general values of r and n follows by the rescaling 

pf---{z) = J^^^vf-^--'^^ {n^r\) , , = 0, 1, 2. (4.30) 

This equality follows from the definitions ()4.10p - ()4.13p changing variables in the integrals. Note 
that this rescaling is consistent with the ODE (|4.2p . □ 



^2(z) = 1 e"V3e-f (2^+i)^,^ei-'^'''-"-^^'' ( 1 + O^z-^l^)] , (4.29) 
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Figure 5: Hankel contour (solid) used for integration in ()4.34p and its deformation (dashed). 



Lemma 4.9. Assuming r = n = 1, the following relations between pj, j = 0,1,2, and qj, 
j = 1, 2, 3, hold 

= 7i^e"'^'^''?i(^)' ^ e C\M-, (4.31) 

"^'^ ^ |_i_e'^"^e° /2z^g;'(z), lmz<0, ^ ^ 

= -:73-e"'"^e"'/'(--)'^92(^) ^ G C\M+. (4.33) 



Proof. We will prove (|i:3T]) - (f03]) for i/ G (-1, oo) \ Z. When G Z the equalities still hold 
true by continuity in u. 

We start with the proof of ()4.3ip using the convergent series representations 



OO ^. „ 

.•^q'Kz) =z''^^ / f'^-i+^eSI^"^ dt. 



Recall the contour integral formula for the reciprocal of the Gamma function [T] 

where the integration is over a Hankel contour as shown in Figure [SJ The Cauchy theorem allows 
us to deform this contour of integration into a contour coming from infinity under the angle —cj), 
avoiding the negative real line, and tending to infinity under the angle (p, for it/2 < (f> < 37r/4, 
see Figure [5i Then 

/"OO 2 1 /"OO /• 2 

/ y'^+'^e-^-^^dy = / y^^^ / ^"^-^-'=6-'^-°^+* dt dy, 

Jo 27ri Jo J 



where the inner integral is taken over the deformed Hankel contour. In the inner integral (thus, 
for fixed y > 0) we change variables t = y/s 



. .,2 1 /"OO r . y2 



1 2 1 



r(i/ + A: + 1) Jo 27ri Jq Jrj; 
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Figure 6: Contour of integration F* (solid) used in the proof of Lemma 14.91 and the contour Fi 
(dotted). 



The new contour F^[ of integration emerges from zero under the angle —(j) and ends in zero under 
the angle (p, where tt/2 < cj) < Stt/A. It has counterclockwise orientation, see Figure [6l Using 
the Fubini theorem we can change the order of integration to become 

1 /"OO 2 If f°° 2 

' y^+V^-"^dy = — / s'^-i+M e-^-"^+^dyds 



F(z^ + /c + 1) 27rz 



r* J-oo 
1 

27ri jp* 



2 

e-V-^J'+f dyds. (4.35) 



We complete the square in the first term of the right-hand side of (|4.35p and evaluate the 
Gaussian integral. Then we can deform the contour F* into Fi. Hence the first term of the 
right-hand side of ()4.35p equals 

1 si r .^_i+fc -L 



e-2 / t'^-^+^e^"^ dt. 



The second term in the right-hand side of (j4.35p vanishes, since we can contract the contour of 
integration thanks to Lemma 14.101 This completes the proof of ()4.3ip . 

Next we prove ()4.33p . We claim that the following chain of equalities holds 

P2{z) = poi-z) = -^e^y\-zrq'li-z) = ^e^^^^-'^) e'^' '^-zY 4{z), (4.36) 



for z G C \ M^. Here, a tilde means that the parameter a on which the quantity depends has to 
be replaced by —a. Then (|4.33p is immediate from the claim. We now discuss the equalities in 
(|4.36p . The first one follows from (|4.10p and (|4.13p . The second equality relies on (|4.3ip . The 
last equality is a consequence of ()4.2ip for j = 1, 2, where one has to take special care about the 
position of the branch cut. 

Finally we prove (I02]) . Using ()CTI]) - (fiT3D and ([LS]) we find 

E T^i^ — -(zy)'^e-J^V2-«. 

This function belongs to the one-parameter family of entire solutions of the ODE (j4.2p . Therefore 
it is characterized by its value at zero 

27f „,2 



F(l - v\ 



The idea is now to build an entire function as a linear combination of the 

-z^<(-z), J = 1,2,3,4 

and find its value at zero to obtain the desired equality. 



24 



When x > we have 



2 



x'^iq'iix) + q'iix) - 4{x)) = / s^^^e^-"--^' ds, (4.37) 



where we used (j4.2ip and made the change of variables s = xt. The integral is taken over the 
dashed contour in Figure [5j Clearly the right-hand side of this expression is entire in x. Therefore 
the equality is valid in the right half plane {z G C | Re z > 0}. The analytic continuation to the 
left half plane takes the form 

-z^ie^^'-'q'liz) + q'iiz) - q'.^{z)) = [ s^-^e^-'^'^^' ds, Re z < 0, 



where the integration is over the same deformed Hankel contour and the argument of z is chosen 
in the interval (— 37r/2, — 7r/2). 

We can use the right-hand side of (|4.37p to determine its value at zero 

2m 



where we also used (|4.34|) . Hence we find the equality 

1 ai 

-2is\n{v'i:)pi[z) + Pq{z) + P2{z) = —?=-& ^ z'' [ql{z) ^ q[{z) - q'^{z)), Rez > 0. 



Now apply ()4.3ip . (|4.33p (assuming Im z > 0), and ()4.22p to obtain (|4.32p . In the other quadrants 
similar relations can be found. □ 



Lemma 4.10. The function 



fa{s)= f e-^-"^+f dy 



oo 



is analytic in C \ {0} and satisfies fa{s) = 0{s) as s — )■ uniformly within the sector \ args| < 
3-7r/4 — 5, for any small 6 > 0. 

Proof. We rewrite fais) as 

/„(s) = e^y|erfc(-^), * = ^ " «, 
where the complementary error function is defined as 

erfc(2;) = —= / dy. 
This function has uniform asymptotics 

erfc(z) = ^(l + 0(z-2)), 

as z — )• oo such that | arg z\ < Sir /A — 6, for a small 5 > 0, see Then fa{s) has the asymptotics 
as s — 7- such that | arg s\ < 37r/4 — 5. This proves the lemma. □ 
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4.4 Wronskian matrix 



From the three special solutions pj, j = 1,2,3, introduced in the previous section we construct 
the Wronskian matrix 



fpoiz) Piiz) P2{zf 
Wn{z)=\p'^{z) p[iz) p'2(^) I xdiag(l,l,e±^-*), 
Viiz) p'iiz) p'^{z)^ 



±\m.z > 0. 



(4.38) 



This Wronskian matrix will be used in the first transformation of the RH problem for Y{z), see 
Section 15.11 

Proposition 4.11. The Wronskian matrix Wn satisfies the following RH problem. 

(1) Wn{z) is analytic for z G C \ M. 

(2) Wn has jumps on M+ and given by 



diag ( 1, 
diag 



1 
-1 1 







1 



WnA^)=WnA^) X { 



(3) As z — )• oo, we have that 



Wn{z) = 3"^r'^et^^z^diag(n"V-^/3zl/^l,nr^/3z~^/^) 



X £ 



(4.39) 



'0 * 0> 



^0 *^ 



X /+ *\z-^/^+\* + 0(^-^1 



0> 



,0 * Oy 





a;2 




)(i 


1 











X diag(l, —a ^ , a)e 



ne{z) 



(4.40) 



as z ^ oo in closed sectors of the upper half plane {z £ C \ Imz > 0}. Here the * denote 
certain constant real entries and 

(4.41) 



We have the same expression as z ^ oo in closed sectors of the lower half plane, but then 
with Lo replaced by uj~^ and diag{l,—a^^ , a) by diag{l, a, a~^). 



Proof. This follows from Lemmas 14.51 and 



□ 



Note that Wn{z) is not uniquely determined from the above RH problem, since we did not 
specify the behavior near the origin z = 0. This will be done in ()5.9p below. 

Lemma 4.12. The asymptotics ofWn can be rewritten as 

Wn{z) 



u+2 2v-2 

-ir 3 esnz s 





n 


* 






1 





1 


* 


+ 0{z-') 




^0 










diag(n-V-4/3^^/3^ 1^^^4/3^-1/3) 

X A± diag(l, cjT\ cj±i)e'^®("), (4.42) 



as z —7- oo in closed sectors of the half plane {z G C | iblmz > 0}, with a as in ()4.4ip and 



1 -uj^ ui 
1 1 



A_ :-- 



-U UJ 



V3 U 



1 1 1 



(4.43) 



a; UJ 
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The prefactors in ()4.43p are chosen such that both matrices have determinant 1. 

Remark 4.13. It is easily seen from ()4.2p and ()4.38p that the determinant of Wn satisfies the 
following linear differential equation 

zidetWnYiz) = {2iy -2)detWn{z), zeC\R-. (4.44) 

Hence, we have for a certain nonzero constant K that 

det Wniz) = Kz^''-^. (4.45) 

5 Steepest descent analysis for Y{z) with quadratic potential W: 
regular cases 

In this section we will perform a Deift/Zhou steepest descent analysis for Y{z) in the situation 
that W is quadratic (|2.12p and the triplet {V, W, r) is regular in the sense of Definition 12.51 This 
analysis consists of a series of invertible transformations and results in the proof of Theorem l2.6l 

5.1 First transformation Y X 

The idea behind the first transformation Y X is inspired by [22]. We will multiply the RH 
matrix Y on the right with the inverse transpose of the Wronskian matrix Wn- Recall from 
(|4.38p that Wn is constructed from the solutions to the third-order differential equation satisfied 
by ho,n- 

In the analysis we will work with the alternative system of weight functions in (j2.6p rather 
than those in (|2.5p . Thanks to (|2.7p . (|4.3p . and ()4.10p these weight functions are 

'wo,n{z) = e-^^('^Po{z), 

< wi,n{z) = e-^^('hp'o{z), (5.1) 
W2,n{^) = e-"^(^)(Vo'(z) + (1 - u)p'o{z)). 

Note that 

wi^n <ind W2^n have been interchanged with respect to (j2.7p . The labeling (jS.ip will 
be more convenient for us. Defining 

/I \ 

D{z) = z l-u\, (5.2) 
\0 z ) 

(|5.ip leads to the equality 

(z) W2,n{z)\ 

e-^n^)W^{z)D{z) = i * * * , (5.3) 

\ * * * / 

where * denotes unimportant matrix entries. This identity motivates the transformation Y t-^ X 
given below. 

Definition 5.1. We define X by 

X{z) = z^'^P-^Y{z) diag (l, D-i(z)iy-^(z)e«®(^)) , (5.4) 

for z G C \ M, where Q{z) is defined in ()4.15p and D{z) in ()5.2p . P„ is an invertible matrix to 
be specified below (|5.7p . 
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We claim that X{z) is the unique solution to RH problem 15.21 below. Note that the jump 
on the positive real axis is simplified at the cost of creating a jump on the negative real axis. 

RH problem 5.2. The 4x4 matrix-valued function X : C\M — )• C^^^ defined in ()5.4p satisfies 
the following conditions. 

(1) X{z) is analytic for z G C \ M. 

(2) We have 

X+{x)= X„ (x) Jx{x), xeR, 
where the jump matrices Jx are given by 



Jx{x) = diag ( l,e 



-ne-{x) 



A 


^-nV{x) 





0^ 





1 














1 


1 










1/ 



diag l,e 



ifx>0, 



Jx(x) = diagfl,e-"®-(^) 



























-1 







^ 











diag l,e 



if X < 0. See (j4.15p for the definition of Q. 

(3) We have as z —)■ oo with itlm z > 0, 

X{z) = [I + Oiz-^)] diag (z", z- V , z-^, z-t) diag {z^'\ z-^'^A^^^ 

X diag (1,1, cT=^So-^^) , (5.5) 

where a and A± are defined in (j4.4ip and ()4.43p . respectively. 

(4) z 0, z G C \ M, we /ia?;e 



X(z)diag(|z|-'^/Mz|'^/Mz|-'^/Mz|^/2) = 0(1), i/i/ > 0, 
X(z)diag(l,(log|zri,l,(log|zri) = 0(1), ifv = Q, 
X{z) = 0(\z\'l^), X-^{z) = 0{\z\^l^), if-l<u<{). 



(5.6) 



Proof. The jumps of X follow directly from ()4.39p . ()5.3p and the definitions. Next we check the 
asymptotics for z — t- oo. From (|2.9p and (|5.2p we obtain 

Y{z) diag (1, D"^(z)) = (/ + 0{z-^)) diag (^z", z-"/^^ z-"/^-!^ ^ 
as z — )• cxD. Moreover, Lemma l4. 121 vields 

z-/2w^-^(z)e"e(^) = ir-'^'e-^zi- 

X diag (r^/3nz-i/3, 1, r-^/^n-iz^/^) A^'^diag (l,c7±\a^^) 












0^ 




G 


1 


* 


1 





+ 0(z-) 






V* 


* 
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as z — )• oo in a sector of itlm z > 0. Combining this we get 



z^l^Y[z) diag (1, L>-i (z) 1^-^(2)6"®^^) 



/I 







2-1/ 

ir 3 e 


3n n 


* 







2+1/ 
3 


Vo 




* 




n + l 

Z 3,2; 







* 




, , - u -I 

IT 3 6 3nn~ / 



as z — )• oo with itlmz > 0, for certain constants *. Finahy, we define P„ to be the constant 
matrix in the square bracket in (|5.7|) . which yields (|5.5|) . 

Next we show the behavior near the origin in ()5.6p . Recah the three weight functions 
e-"^(^)/io,n(x), e-'^^(^)x/i'o „(x), and e-'^^(^)(x/i^'_„(x) + (1 - v)h!^ Jx)) in the RH problem for 
Y . It follows from the explicit form of Y that 



y(z) = o(i), ifi/>o, 

y(z)diag(l,(log|z|)-i,l,(log|z|)-i) = 0(1), if z. = 0, 

diag(l, Iz]-'^, Iz]-'^, |z|-^) = 0(1), if -1< < 0, 

y(z)-^diag(|z|-^ 1,1,1) = 0(1), if -1< < 0, 



(5.8) 



as z — )• 0. On the other hand, we see from the definitions of pi, i = 0, 1, 2, that there exist some 
constants A, i?, C, D that depend on v such that 



po(^) = ^^"+o(kr+'), 

f5z'^ + 0(|z|'^+i) + C + 0(z), ifz^^O, 

|5 + o(z) + iogz((:; + o(z)), ifz^ = o, 



pi{z 



(5.9) 



as z — 7- 0. This, together with ()5.2p and ()4.38p gives 



^z^' + C'dz 



Odzl'^) + 0{l) + 0(log |z|) Dz'^ + 0{\z\''+^) 



D^{z)Wn{z) = ( uAz" + 0{\z\''+^) 

0{\zn 



Oi\zn + Oiz) 
0(|zr) + 0(l) + 0(log|z| 



i^Dz^ + Oilz 
0{\zn 



as z — 7- 0, where we understand that the ©(logjzj) terms are absent if 7^ 0. Also in case 
1^ = the (1, 2) and the (3, 2) entry have to be replaced by 0{z). Note that the terms of order 
Odzl*^"^) in the last row all cancel. 

Using the cofactor formula for the inverse transpose and the fact that 



det(L»^ (z)VF'„(z)) = Kz^", 



see (j4.45p . we find 

D-\z)W-^{z) 
1 



0{\z\^^) + 0{\z\^) 0{\z\^^) 
0(|z|2^) + 0(|z|'^) + 0(log|z|) Odzp'^) 
0{\z\^^) + 0{\z\') 0{\z\^^) 



0(|z|2-) + 0(|zr) 
0(|z|2^) + 0(|zr) + 0(log|z|) 
0(|z|2-) + 0(|zr) 



where again the ©(log |z|) terms are absent if 7^ 0. Using this with (|5.8p and ()5.4p . we arrive 
at ()5.6p after a straightforward calculation. □ 



29 



5.2 Second transformation X ^ U 

The second transformation X ^ U serves to (partly) normalize the behavior at infinity. To 
do this we will use certain functions related to the vector equilibrium problem. We call these 
functions A- functions and define them as a transformed version of the A-functions in [23] . More 
precisely, denoting with Aj^'~^(z), j = 1, 2, 3, 4, the A-functions used in ^23j, we will work with the 
square root versions \j{z) = 2)^^{^/z). 

Lemma 5.3. There exist functions Xj, j = 1, 2, 3, 4, analytic on C \ M that satisfy the following 
conditions. 

(a) As z ^ oo we have 

Xi{z) = V{z)-log{z)-£i + 0{z-') , 

X2iz) = ei{z) + \ log(z) + + 2?Z-2/3 ^ Q (^-1) ^ (5 

C \ \og{z) + Cuz-'^/'^ + Duj'^z~^l^ + O (z-i) Im z > 0, 



Xziz) = e2{z) + 



\ log(z) + Cuj'^z-^/^ + Duz-"^/^ + O {z'^) Im z < 0, 
\ \og{z) + Cw^z-i/s + 1)^2-2/3 + (^-1) im z > 0, 



A4(Z) = 6*3(2) +< Y . . 1 /q O O/Q / IN ' (5-11) 

^' \\\og{z) + Cujz-^l^ + Duj^z-^/^ + 0{z-^) lmz<0, 
where C, D and ii are real constants. 
(b) There exists a positive integer N , two sets of ordered numbers 

= bo < ai < bi < 02 < ■ ■ ■ < On < bN < cat+i = 00, 

1 = ao > ai > ■ ■ ■ > qat = 0, 

and constants C2 > and C3 > such that the X-functions satisfy the following jump 
conditions: 

(i) On we have 

Ai,±-A2,^ = on {aj,bj), j = 1,...,X, (5.12) 

Ai,+ — Ai^_ = —2TTiaj on {bj,aj+i), j = 0, . . . , X , 

A2,+ - A2,- = 2iTiaj on {bj,aj+i), j = 0,...,X, 

A3,±-A4,=F = on (ca, 00), (5.13) 

A3,+ - As - = -27ri/3 on (0,03), 

A4^+ — A4^_ = 2iTi/3 on (0,03). 

(ii) On R~ we have 

" ■ onM", 
A2,± - A3,=p = ±27ri/3 on (-00, -C2), (5-14) 



'1,+ 


-Ai,- 


= -2TTi 


'2,± 


-A3,T 


= ±27ri/3 


'2,+ 


-A2,_ 


= 27rV3 


'3,+ 


-A3,- 


= 27ri/3 


'4,+ 


-A4,- 


= 27rV3 



A2,+ - A2 - = 27ri/3 on.(-C2,0) 
A3,+ - A3 - = 27ri/3 on(-C2,0) 
Aa + — Azi _ = 27ri/3 onM^. 



Proof. Define 



Xj{z) = 2Xf^{^/z), j = 1,2,3,4:, zeC\M, (5.15) 
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where we take the potentials as in ()2.19p . Then (a) fohows from [23] Lemma 4.14] and ()4.16p . 
Note however that error terms stated there are not optimal. (b)(i) is direct from j23j, Lemma 
4.12]. To prove (ii) we need some preparations. By \23\ Definition 4.1], it follows that for x > 

gf''iix)-gf''{-ix) = 7Ti, 
g2,±iix)-gl'^i-ix) = 27Ti/3, 
gf''iix)-gf''i-ix) = 7ri/3, 



and 



5^ (-^) = (^) T 2^V3, ± Im z > 0, 

gf'^i-z) = g^'^iz) T vri/3, ± Imz > 0. 

Then [23, Definition 4.11] and (liTTll yield for x > 



X^'^iix) 


- (- 


-ix) 


= —iri, 




X^^_{^x) 




-ix) 


= ±Tri/3, 


C^^ <X < +00 




-A^^(- 


-ix) 


= 7ri/3, 


<x < c^^, 


X^^_{^x) 




-ix) 


= 7ri/3, 


<x < c^^, 




- Ar(- 


-ix) 


= 7ri/3. 





This, together with (|5.15p . implies (ii) in (b), where we take C2 = {c2^)'^. □ 

The A-functions also satisfy a number of inequalities stated in the following lemma. 
Lemma 5.4. 

Re(A2,+ - Ai,_) < 0, on{bj,aj+i), j = 0,...,N, (5.16) 

Re(A2,+ - As,-) < 0, on{-C2,0], (5.17) 

Re(A4,+ -A3,-) <0, on[0,C3). (5.18) 

Proof. This follows from ()5.15p and \23\ Lemma 4.13]. □ 

Remark 5.5. For the constants C2 and C3 in Lemma |5.3[ we have C2 > —y*{a) and C3 = if 
a > 0, and C3 < x*{a) if a < 0, where x*{a) and y*{a) are given in Lemma 14.61 The functions 
Ai and A2 are defined and analytic on C \ (—00, ^at], whereas A3 and A4 are defined and analytic 
in C \ M. 

For future reference we specify the behavior of the A-functions near the origin. 

Lemma 5.6. In a neighborhood of the origin the X-functions, defined in Lemma \5.3\ have the 
following behavior for 2; — 0. 

(a) In Cases I and IV there exists a constant ci > such that 

(Ai - X2){z) = =F2vri ± iciz^/^ + 0{z), ±lmz > 0. 

(b) In Case III there exists a constant C2 > such that 

iX2 - X3)iz) = ±^ + C2z'^/^ + 0{z), ±lmz>0. 

(c) In Cases I and II there exists a constant C3 > such that 

(A3 - Xi){z) = ± ic3Z^^^ + 0{z), ±Imz > 0. 
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Proof, (a), (b), and (c) follow from ()5.15p and the proofs of ]23\ Lemma 7.5, Lemma 7.1, Lemma 
7.2] respectively. □ 



We can now define the transformation X i— )• [/ in a similar way as in 
Definition 5.7. We define the 4x4 matrix-valued function U by 

U{z) = [/oe"^X(z)e"^(^), (5.19) 
where Uq is a constant invertible matrix to be determined later, 

G = diag (Ai - A2 - 01, A3 - 02, A4 - 63) , 
and L = diag(^i, 0, 0, 0) with £1 the constant from Lemma [ 



RH problem 5.8. The matrix-valued function U defined in (|5.19|) is the unique solution of the 
following RH problem. 

(1) U (z) is analytic for z £ C\M. 

(2) We have 

Tj (^\ _ TT j diag ((Jc/)i (x), {Ju)3 (x)) , X G E+, 

u+[X) ^-W ^ cliag(e^"S(J[/)2(x),e-'^"% xGR_, 



where 



g«(A2,+ -A2,-) J ' 

g 'i/7rign(A2,+ -A2,_) q 

_g"{A2,+ -A3,_) e'^'^^e 
pn(A3,+-A3,_) n(A4,+-A3,_) 

iJu)s • 



{Ju)2 - \ „n{A2,+ -A3,-) pUiTipn{\-i,+^\-i^-) I ' 



gn(A4,+-A4-) J ■ 

(3) We can choose the constant matrix Uq such that 

U{z) =[l + {z-^)] diag (l, z^l'^ .z'^l'^ , 1^ 

X diag [z^l'^, z-^'I^A^^ diag (l, 1, ,&^^) , (5.20) 

OS z —7- 00 with ibim 2 > 0, where a and A± are defined in (|4.41|) and (|4.43|) respectively. 

(4) U{z) has the same behavior as X{z) near the origin; see (|5.6|) . 

Proof Noting that e"(^i'+"-^i'-) = e"(^4'+~^'' -) = 1 on M~, due to Lemma l^^ b) and our as- 
sumption that n = mod 3, the jump condition for U in item (2) follows from a straightforward 
calculation. 

To show the asymptotic behavior of U for large z in item (3), we see from Lemma l5.3f a) 
that 



gnG(^) ^ j-j^g (^^-ng-n^l ^ J + ^ \ Z 





J^±, 



+nE + diag (O(z-i), 0(^-1), 0(^-1), ©(z-^)) 



as z —7- 00 with iblmz > 0. Here E = D + C'^/2 and 

ri+ = diag(l,a;,a;^), ri_ = diag(l, w^, w). (5-21) 
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Then, by the asymptotic behavior of X in ()5.5p and the definition of U in ()5.19p . we get 
U{z) = Uq{I + diag (l, zV3^ ^-1/3^ 1^ 

X diag [z''/^,z'''/^Af'^ 



I + nC 




17± 



'^/'^ +nE 









-2/3 



+ diag(0(z-i),0(z-i),0(z-i),0(z-i))] xdiag(l,l,a±\a^i), 

as z — oo with itlmz > 0. We can move the terms within square brackets to the front at the 
expense of an extra constant contribution. This follows from the observation that 



^0 1^ 



^0 1 0^ 



A^' 9.±A\ =110 1 and A;^' ^t^± = I 1 I • 



(5.22) 



.0 1 0, 



.1 0; 



This gives 

U{z) = Uo 





A 















1 


* 


* 










1 







Vo 





* 


1/ 



diag (l, z^/^, l) diag (z'^/^^ z-'/^A^^ 



X diag(l,l,(j=^\cj'Fi), (5.23) 



as z — )• oo with itlmz > 0. By setting Uq to be the inverse of the constant matrix between the 



above square brackets, we arrive at the claim in item (3). 

Recalling our assumption that n = mod 3, we find from Lemma 15.3( b) that 

jn(Ai, + -Ai,_) 2 \ 



□ 



iJu)i 



g"-(A2,+ -A2,-) 

^-2mnaj gn(A2,+-Ai,_)\ 
^2mnaj j ' 



L, . . . , X T , 



, on {aj,bj), j = 1, 

on {bj,aj+i), j = 0,...,N. 



(5.24) 



(Ju) 



-i^7rign(A2,+ — A2,-) 
-1 



gj^7rjgra(A3,+-A3,_) ) ' 



U)3 



_g'T.(A2,+-A3,_) ^uni j ' 
'l gn{A4,+ -A3,_)\ 

,0 1 j ' 

gn(A3,+-A3,_) ^ 



on (-00, -C2), 



on (-C2,0), 



(5.25) 



on (0,03) 



on (,C3,ooj 



(5.26) 



e"(-^4,+-A4__) ) ' 

Note that the diagonal terms of {Ju)i are 1 on (0,ai) and {bN,oo) 



5.3 Third transformation U t-^ S 

In the third transformation we open lenses around [jf=i{aj, bj), (—00, — C2) and (ca, +00), which 
are denoted by Li, L2, and L3, respectively. See Figure [7] for a plot of the lenses in Case IV if 
= 3. These lenses are chosen such that the following estimates hold. 

Lemma 5.9. We can find a neighborhood Li of \J^^^{aj,hj) such that 

Re(Ai - A2)(z) < /or zGLi\M+. 
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Figure 7: Jump contour for S consisting of the lips of the lenses Li, L2, L3, and the axes. 
This is a picture of a Case IV situation with = 3 since ai = and 02,03 > 0. The lips of the 
lens around supp(/i3) intersect the lips of the lens around (a^, 5^) if and only if C3 E {a^, 



We can find a neighborhood L2 of (— cxd, — C2) such that 

Re(A3 - A2)(z) < forz£L2\Mr, 

and such that 

{z e C I Rez < -i?, |Im(2;)| < -5Re(z)} C L2 

/or some e > and R > 0. 

Finally, we can find a neighborhood L3 of (03, +00) such that 

Re(A3 - A4)(z) < forz£L3\R+, 

and such that 

{z £ C \ Re z > R, |Im(2;)| < eRe (z)} C L3 

for some e > and R > 0. 

We can assume that Li, L2, L3 are symmetric with respect to complex conjugation. Moreover, 
we can assume that the intersection of the boundary dLj with the upper (or lower) half plane is 
a collection of smooth curves, which can be oriented so that the real part strictly increases. We 
will often refer to dLj as the lips of the lens Lj . 

Proof. In view of the definition of A- functions in ()5.15p . the statements about Li, L2 and L3 are 
immediate from Lemmas 7.1, 7.2, and 7.5 in i23]. □ 



To define the third transformation, we observe that in ()5.24p - (|5.26p we can factorize 
{Ju)2=[^ 1 Jl-1 oJio 1 J' (5.27) 

^Juh = (gn(A3-A4)- 1) J) (en(A3-A4)+ ' (c3,Oo), 

with the aid of (j5.12p . (j5.13p . (j5.14p and the fact that n is a multiple of three. 
We then define the third transformation [/ 1— )• S by successively setting 

r /±e±''^*e"(^3-^2)(^)^2,3, forzGL2nC±, 

T{z) = U{z) X <^ / ^ e"(^3-A4)W^^^3^ ' for z G L3 n C^, (5.28) 

[ /, elsewhere, 

and 

SU) = / X T e"(^i-^2)(^)i?2,i) , for z G Li n C±, 

\ T{z), elsewhere, 

where denotes the upper or lower half plane respectively, and Eij denotes the 4x4 elementary 
matrix of which all entries are 0, expect for the (i,j)-th entry, which is 1. Then S is the unique 
solution of the following RH problem. 
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RH problem 5.10. The matrix-valued function S satisfies the following conditions. 



(1) S{z) is analytic for z £ C \ S5, where is the contour consisting of the real axis and the 
lips of the lenses Li, i = 1,2,3, see Figure^ 

(2) For z E S5, S has a jump 

q ( \ =Q m/ ^^^S (('^■5)1 (^)' (•^•5)3 (^)) , for z inR+ and the lips of Li,L3, 
- :y-{z) I ^.^^ j^^^i^ ^j^^^ {z),e-''^') , for z in M" and the lips of L2, 



where 



(J, 



s)i 



^—27vinaj gn(A2,+ — Ai,_) 
gn(Ai-A2) Y \ ' 



and 



_g"(A2,+-A3,_) gi^vri ^ ' 

A 
-1 Q) ' 

1 



2 gn(A4,+-A3, 

1 

r 

-1 

1 0\ 

gn(A3-A4) ^ J ' 



on {aj,bj), j = 1,...,N, 

, on (6j,aj+i), i = 0, . . . ,iV, 

on the lips of Li, 
on the lips of L3, 

on (-C2,0), 
on (-00, -C2), 

on the upper /lower lip of L2, 

on (0,03), 

on (c3,+oo), 

on the lips of Li, 

on the upper/lower lip of L3 . 



(3) S{z) has the same behavior for 2; — t- 00 as U{z). 

(4) S{z) has the same behavior near the origin as U{z) (and X(z)), see (|5.6p . provided that 
2; —7- outside the lenses that end in 0. 

Proof. The jump condition is straightforward by ()5.27p and item (2) in RH problem 15.81 The 
asymptotic behavior of S follows from the definition in ()5.28p - (|5.29p and the large z behavior 
of the A-functions in (ISTOjl - ffSlTTl . □ 

In view of Lemmas 15.41 and 15.91 it is easily seen that each entry of Js is either constant or 
exponentially small as n — )• +00. 



5.4 Global parametrix 

In this section we look for a global parametrix S^°°h This will be a good global approximation 
of the matrix-valued function S when n is large. Ignoring all exponentially small entries for 
n — )• 00 in Js, we are led to the following model RH problem for S^°°\ 

RH problem 5.11. We look for a 4x 4 matrix-valued function S'(°°) that satisfies the following 
conditions. 
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(1) S^°°\z) is analytic for z G 

(2) For xeR, 3^°°^ has a jump 



m 



diag (e*"^*, (J5(oo))2 {x),e , /or x in . 



where 



' f r 

diag (e-2™"j , e2™a,^ ^ (6^,0^+1), j = 1, . . . , - 1, 



on {aj,bj), j = 1,...,N, 



on (0, oi) U {bN, +00), 



and 



diag (e-'^^\ 6^^^% on (-C2,0), 

OA , . 

on (-00, -C2), 



on (0,03), 



°°) J3 



-1 0; ' 

r /2, 

1 

-1 



, on (c3, +00). 



(3) 2; — )• 00 iblmz > 0, we have 

S^^\z) = [l + [z-^)] diag (1, 1) diag z-''/*'^^^ 

X diag(l,l,(7=^\o-=Fi) . 

In this section we will construct a solution to the above RH problem. 

Transforming the global parametrix: 5'^°^^ 1— )■ 

We will look for a solution S'(°°) to RH problem 15. Ill in the form 

= iV^(z) diag (^'^,2^/3) diag (1,1, cj^\cj^^), ±lmz>0. (5.30) 



Then Ny must satisfy the following RH problem. 

RH problem 5.12. We look for a 4 x 4 matrix-valued function Njj satisfying the following 
conditions. 

(1) Ny{z) is analytic for z G C \ M. 

(2) For X G M, Ny has a jump 

Ny,+{x) = Ny,-{x) 



diag ((JatJi (x), ( JArJg (x)) , for x in 
diag (-1, {Jn,)2 (x), -1) , for x in 



where 



2v-2 ' 

z— 

■2-2v 

-z^- 



on {aj,bj), j = 1,...,N, 



diag (e-2™°. , e2™".) , on (6^,0^+1), j = 1, . . . , iV - 1, 
h, on (0,ai) U (6Ar, +00), 
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Figure 8: Plot of the Riemann surface TZ for Case V (i.e. ai, C2, C3 > 0) and genus g = 2. 



' diag (^e-t(i+2.)i^g|(i+2.)i^ {-C2M, 

OA , , 

on (-00, -C2j, 



and 



('^AfJ3 = 

(3) As z — 7- 00 mt/i ± Imz > 0, we have 



diag 

1 
-1 



1 0, 

g-f(l+2i.)i 



-63 



f (l+2z/)j 



on (0,03), 
on (c3, +00). 



N^{z) = + diag 



M2 ^1/6 ^-1/2 



,z-i/6)diag(l,^±^). (5.31) 



Boundary value problem on a Riemann surface 



The next step in the construction of the global parametrix is to find a 'Szego function' on a 
certain Riemann surface. We define a four-sheeted Riemann surface Tt as follows. We let T^j, 
j = 1, 2, 3, 4, denote 



AT 



7^1 



k=\ 



7^2 = c\ ^|J["'^''^^]^(-°°'-^2]j 



7^3 = C \ ((-00, -C2] U [C3, 00)) , 



7^4 



C\ [C3,00). 



We connect the sheets T^j, j = 1,2,3,4, to each other in the usual crosswise manner, e.g. TZx is 
connected to 7^2 along the cuts [a^, 6^], A; = 1, . . . , A^. The Riemann surface is compactified by 
adding two points at infinity: ooi is added to the first sheet while 002 is common to the other 
sheets. We define B as the union of four small disks, one around the origin of each sheet. We 
denote by C C 7^ the contour consisting of the intervals [ai,6Ar] on the first two sheets, [— C2,0] 
on sheets 2 and 3, and [0, C3] on sheets 3 and 4. See Figure [8] for an illustration in Case V. 

We want to construct a Szego function on this Riemann surface, i.e. we look for a scalar- 
valued function f onTZ satisfying a boundary value problem. We denote the restriction of / to 
the jth sheet with fj{z), j = 1,2,3,4. 

Boundary value problem 5.13. We look for a scalar-valued function f satisfying the following 
conditions. 

(a) f is analytic onTZ\C, and there exist constants Cx,C2 such that < Ci < |/(^;)| < C2 < 00 
onn\{CuB). 
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(b) f has the following jumps on C : 





2v-2 


X G 


iaj,bj), 


(5.32) 






X G 




(5.33) 


/2,H-(X) 




X G 


{bj,aj+i), 


(5.34) 






X G 


(-C2,0), 








X G 


(-C2,0), 




/2,±(X) 


= h,Tix), 


X G 


(-00, -C2), 




/3,+ (x) 




X G 


(0,C3), 






= /4,-(x)e-f(i+2'^), 


X G 


(0,C3), 




/3,±(^) 


= At (2^)' 


X G 


(C3,00), 





where in (|5.33p - (|5.34p we have certain real numbers 13 j, j = 1, . . . , N — 1, to be specified 
in Lemma \5.17\ 



(c) f is regular at z = 00 in the sense that 

h{z) = c + 0{z-^] 

and 





-1/3 



UJ I + C32; 



-2/3 




+ 0{z 



ci /O, 



as 2; — 7- 00 in the upper half plane. 



(d) f has the following behavior around the origin of each sheet for Imz > 0. The behavior 
depends on the particular case we deal with, see Section\2.4\ 



/i(^) 



f /^oe'?(^-''')^3(2--i) (1^0 (^1/2)) ^ Cases I and IV, 
\o{l), in Cases II, III, and V, 



Koe-x(i-2-)^i^(5-2.) (1 + 0(^1/2)), 

KlZ^(l+2^)(l + 0(z)), 
^K2Z-T2i^+^-) (1 + 0(^1/2)), 

{K^e-'r^^^+^^h-U^+^-) (1 + (zV2)) 
K2ef(i+2.)^-J,(i+2.) (1 + 0(^1/2))^ 

^K4ef(i+2-)^-i(i+2-)(i + o(^)), 

/^3eS(l+2-)z-li(l+2-) (1 + O (zl/2)) , 

^5e-f(i+2'^)zi(i+2-)(l + 0(z)), 



in Cases I and IV, 
in Cases II and V, 
in Case III; 

in Cases I and II, 

in Case III, 

in Cases IV and V, 

in Cases I and II, 

in Cases III, IV, and V, 



(5.35) 
(5.36) 

(5.37) 
(5.38) 



as z — 7- with lm.z > 0. The behavior in the lower half plane can be obtained using the 
symmetry condition f{z) = f{z). The kj, j = 0, . . . ,5 are real constants. 

Note that (c) is not an extra restriction as it is implied by (a). 

We will settle the solvability of this boundary value problem at the end of this section. 
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Transforming the global parametrix: Ny i— )• My 



Assuming the solvability of the boundary value problem for the Szego function / we can further 
reduce the RH problem for N^, i.e. we look for a solution to RH problem 15.121 in the form 



Ny{z) = CM^(z)diag 



1 



1 



1 



1 



h{zy h{zy HzY u{z) J ' 



(5.39) 



where C denotes the explicit constant matrix 



C 



/c 0\ 

Ci C3 C2 

ci 

VO C2 ci/ 



(5.40) 



with c, ci,C2,C3 the constants from condition (c) in Boundary Value Problem 15.131 The matrix 
C does not influence the jumps but will serve to get the appropriate asymptotics of My{z) for 
large z. Putting 5j := otj + Pj/n for j = l,...,iV — 1, the matrix- valued function My must 
satisfy the following RH problem. 

RH problem 5.14. We look for a 4 x 4 matrix-valued function My satisfying the following 
conditions. 

(1) My{z) is analytic for z G C \ M. 

(2) For X GM, My has a jump 

(My), (X) = {My)_ (X) 



diag ((JmJi (x), {Jm^)^ (x)) , for x in 
diag (-1, (JmJ2 (3^)' -1) ' ^ 



(5.41) 



where 



A 
-1 

diag(e~ g27rmaj' 



and 



1 

-1 

-h, 

1 
-1 



on iaj,bj), j = 1,...,N, 

on {bj,aj+i), j = 1,...,N -1, 
on (0, oi) U (^AT, +00), 

on (-C2,0), 
, on (-00,-02); 

on (0,03), 
, on (c3, +00). 



(3) As z — )• cxo and iblmz > 0, we have 

My{z) = [l + {z-^)] diag (z^'\ z^l\z-^l\z~^l^\ diag (l, . 



(5.42) 



Indeed, the conditions (1) and (2) in the RH problem for My are immediate from ()5.39p . RH 
problem 15.121 and Boundary Value Problem 15.131 For condition (3) one also uses the identities 
(IL2I])-(I522]) and ([QO]) . 
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Constructing the global parametrix 

We will immediately prove the solvability of RH problem 15.141 which then in combination with 
(|5.30p and ()5.39p finishes the construction of the global parametrix. 

We solve RH problem l5.14l bv reducing it to the RH problem for the global parametrix in the 
non-chiral two-matrix model described in [23\ (8.1)-(8.6)]. The latter RH problem was solved 
in the same paper. We denote that solution here as M^'^{z). Note that M^^(z) depends on 
certain parameters n G Z and < < • • • < ctyv-i that will be specified later. We will also 
need the symmetry relation 

M^^(-z) = diag(l,-l,l,-l)M^C(z)diag(l, -1,1,-1), (5.43) 

which is not hard to verify. Now we claim that, 

M^(z) := i^Ldiag(-z^/2j,z^^/^l)M^^(z^/2)diag(-l,-l,^l,±l), ±Imz > 0, (5.44) 

solves RH problem 15.141 Here K is a constant matrix of the form 



K 



/l * * *\ 

1*0 

10 

yo * 1/ 



(5.45) 



for certain constants *, serving to get the correct asymptotics for z — t- oo. Note that an apparent 
2^/^ contribution for z — ?■ oo vanishes due to ()5.43p . L is a matrix of the form 

L = I + Kz-^Es^i, (5.46) 

for a suitable constant k. The matrix L serves to get the correct behavior as z — t- and 
will be constructed in the proof of Lemma 15.151 A straightforward verification shows that 
Ml, defined in (j5.44p indeed solves RH problem 15.141 for the right choice of parameters n and 
< afC < ... < 

Since the RH problem for M^*^ is solvable for any choice of these 
parameters, see [231 Section 8], we have proved the solvability of the RH problem for M^. 

Apart from the solvability of Boundary Value Problem 15.131 we have now finished the con- 
struction of the global parametrix by (|5.30p . (j5.39p . and (|5.44p . We will settle the solvability 
of Boundary Value Problem 15.131 at the very end of this section but first we will discuss the 
behavior of the global parametrix around the origin. 

Behavior of S^°°^ near the origin 

In the next lemma we discuss the behavior of S^°°'^ near the origin. 

Lemma 5.15. The constant matrix L in ()5.46p can be chosen such that S^°°^ has the following 
behavior near the origin z = 0: 

5(°^)(z)diag(zV4,2l/4, ^1/4^^1/4) ^ 0(1)^ (o^^g 

5(~)(z)diag(z-^/^z(l+'^)/^zl/^zV4) = 0(1), m Case II, 

5(°^)(z)diag(z~^/2^z^/^zi/^,z(i+'^)/2) = 0(1), in Case III, (5.47) 

5(°°)(z)diag(zl/4,2l/4^^-^/2^^(l+i.)/2) ^ 0(1)^ (o^gg jy^ 

5(°°)(z)diag(z-^/2,z(i+'^)/2,z~^/2,z(i+^)/2) = 0(1), in Case V. 
Proof. By combining the above transformations (|5.30p . (|5.39p and (|5.44p we get 

5(~)(z) = C7i^Ldiag(-z^/2^ 1, z^'^/\ l)M^^{z'^/^) diag(-l, -1, ^1, ±1) 

xdiag(/fi(z),/2-i(z),a±V3~'(^),^^V4~'(^))diag(z^,zV/3), (5.48) 
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Figure 9: Plot of the Riemann surface TZ and the canonical homology basis 
{Ai, . . . , Ag, Bi, . . . , Bg) for Case V (i.e. ai, 02,03 > 0) and genus g = 2. 



for lb Imz > 0. 

From the construction in [23^ Sec. 8] and the symmetry (|5.43p it follows that 



as z 



in the first quadrant of C, where 



A + Oiz) 



(5.49) 



A 





fa 


a 


* * 




b 


-b 


* * 










* * 




\d 


-d 


* * 




fa 


* * 









* * 


* 




c 


* * 


* 




Vo 


* * 





Cases I, IV, 



Cases II, III, V, 



(5.50) 



where o, b, c, d are certain constants with abod 7^ 0. We then define the matrix L as in (j5.46p with 
the constant k given by k = o/a. The lemma then follows from a straightforward calculation 
using (|5.48p - (|5.50p and the behavior of the Szego functions fi near the origin in (|5.35p -- (|5.38p . □ 

Remark 5.16. By using a finer analysis of the structure in (j5.50p . one can show that each of the 
terms z^^'^'^'^^'^ in the Cases II-V in ()5.47p can be replaced by z*^/^. 



Solvability of Boundary Value Problem 15.131 

It remains to prove the existence of the Szego function / as a solution of Boundary Value Problem 
15.131 On the Riemann surface TZ we construct a canonical homology basis consisting of closed 
curves Ai, . . . , Ag and Bi, . . . , Bg with g the genus of TZ. Here Bj is a closed curve on the first 
sheet going counterclockwise around the union of cuts Ui=i['^fc5 j = ^, ■ ■ ■ , 9- On the other 
hand, Aj is a closed curve on the first two sheets whose intersection with the first sheet is a path 
connecting a point of the cut {aj,bj) to a point of the cut (aj+i,6j+i), and whose intersection 
with the second sheet is the complex conjugate of this path, with the reverse orientation. See 
Figure [9] for an illustration in Case V. 
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Let dwi{z), . . . , dwg{z) be a basis for the space of holomorphic differentials on TZ normalized 
with respect to the cycles Ai, . . . , Ag, i.e., such that 



dwkit) = 6jk, j,k = l,...,g, 

where 5jk denotes the Kronecker delta. 

Let P be a point on the Riemann surface TZ which lies in the upper half plane of one of the 
sheets TZj, j = 1, . . . , 4. There exists a sufficiently small neighborhood U of P, lying entirely in 
the upper half plane of TZj, such that the holomorphic differential dwk allows the representation 

dwk{z) = pk{z) dz, z gU, (5.51) 

for a certain analytic function p^, where we use z as the complex coordinate on the sheet TZj. 
Letting U be the set of complex conjugate points of U lying on the same sheet TZj, we then have 
from the symmetry under complex conjugation that 



dwk{z) = pk{z) dz, z£U, (5.52) 

where we again use z G [/ as the complex coordinate on the sheet IZj. 

Now we construct the Szego function / satisfying (a)-(d) above following Zverovich [36]. In 
the language of Page 135] we must find a solution to Riemann's homogeneous problem where 
^{p) is our function f{z), L is our contour C, and with the divisors D and J prescribing the 
singularities given by 2? = 1 and (using additive notation) 

(mm{^^,^}0i,2 + ^^03,^, Case I, 

+ ^^03,4, Case II, 

=^02,3 + ^04, Case III, 

mm{=^, ^}0i,2 + ^^03 + i^04, Case IV, 

i^Oa + ^^Os + ^04, Case V. 

Thus the only singularities or poles of / are allowed at the points lying over the origin. The 
precise form of the divisor J is due to (|5.35|) - (|5.38|) . 

Denote by G{t) the multiplicative factors appearing in the jump conditions in part (b) of 
the above boundary value problem for /. Note that G(t) satisfies a Holder condition on each 
analytic arc of the contour C, except possibly at the origin of the first two sheets, in the case 
where ai = 0. Thus the behavior at the origin needs to be analyzed separately. 

Following [ini Page 137] we introduce the piecewise analytic function X{q), q £ TZ, as 

X{q) := e^fc^'^G(r)dd^,,,,iT)_ (5_53) 

Here dLUgqg{t) is the discontinuous analogue to the Cauchy kernel, qo C. We will show that 
X{q) solves Boundary Value Problem 15.131 for well-chosen values of the constants j = 1, . . . ,g. 

Lemma 5.17. There exist j3j G R, j = 1, . . . , such that the function X{q) defined in (|5.53p 
solves Boundary Value Problem \5.13[ 

Proof. We first prove that X{q) satisfies condition (b) of the boundary value problem. The 
function X{q) is analytic on 7^ \ (C U Ufc=i ^k) and has the following jumps, see 



X+{t) = X_{t)G{t), foiteC, (5.54) 

X+{t)=X^{t)e~-fc^^G{r)dwUr)^ {oTtGAk, k = l,...,g. (5.55) 
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Note that the jump on A^, k = 1, . . . ,g, is constant. We claim that we can define the real 
constants /Sj such that these jumps are actually trivial. To see this we first show that 



lnGiT)dwkiT), k = l,...,g, (5.56) 

c 

is real. This follows from the symmetry under complex conjugation and the particular form of 
the jump factors G{z). Indeed, using ()5.32p and ()5.5ip - ()5.52p the contribution of the cut [oj, bj] 
on (|5.56p can be written as 

1 /" 21/ - 2 , . , , , , If 2z/ - 2 , , , , . , 
log(x+) dwk{x+) - — I — - — log(x-) dwk{x-) 

3 

1 2u-2 /"^^ 



27r« ./[a,,b,] 3 27ri 



2TTi 



log{x){pk{x) - pk{x)) dx G 



where the first (second) integral is over the interval [oj, 6j] in the upper (lower) half plane of the 
first sheet. Also all contributions 



77^ / \ogG{x)dwk{x), 



are real, where K is any of the intervals (bj,aj), j = 1, . . . , N — 1, on the first or second sheet, 
(— C2,0) on the second or third sheet, or (0,03) on the third or fourth sheet. This follows from 
the fact that on these intervals the meromorphic differential dwk is real and the particular form 
of logG(3;) on these intervals. Moreover, in (|5.56p . the contribution from the gaps (6j,aj+i) on 
the first two sheets is given by (use ()5.33p - ()5.34p ) 



1 

27ri 



{-2TTi/3j) dwk{t) = -(3j5k, 



Hence, the constants f3j, j = 1, . . . ,g, in ()5.33p - ()5.34p can indeed be chosen so that ()5.56p is 
zero. Therefore the function X{q) already satisfies condition (b) of Boundary Value Problem 

Em 

Next we prove condition (d). It follows from the precise definition of G{z) that the quantities 
Xfe) X ill [46 1 Page 138] are all zero except at the points lying over the origin. Hence the function 
X{q) can only have singularities at these points. The order of these singularities is described 
by the divisor £ in Page 138]. In our case this divisor is precisely the inverse of J above: 
£ = . Indeed, this follows from the formulas in [l6] except for the point 0i^2 hi Cases I and 
IV, since there the jump matrices do not satisfy the boundedness and/or Holder conditions. In 
Cases I and IV we have that ai = and C2 > 0, so in the union of disks B the Szego function / 
has the jumps 

/i,±(2;) = /2,t(^)^^, a;>0, (5.57) 

/2,+(x) = /2,-(x)ef (^+2.)^ ^ ^ Q_ (5_58) 

We define the conformal map 

where the cut of the fractional power is chosen along the positive real line, i.e. < argz < 27r. 
This map sends (T^i U 7^2) to a disk B* centered at the origin of the complex ^-plane. Lifting 
the Szego function to the ("-plane we get 



/(C) 



/i(C2) ImC>0, 
hiC) ImC<0. 
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Then ([5371) - ([538]) translate into 



, 4u-4 



/+(C) = /-(C)IC| — , CGMnS*, 
/+(C) = /-(C)e'?(^+'^\ CeiR-nB*, 

where the part of the real line is oriented from left to right and the part of the imaginary axis 
from bottom to top. Applying log at both sides of the equation leads to an additive problem 
that can be solved using Cauchy transforms and the Sokhotski-Plemelj formula 

log/(C) = ^^/_'^^dt + ^^y"°^dt + MC), C(^B*\{Mum~), (5.59) 

where h is meromorphic in B* with no poles except possibly at zero. The asymptotic behavior 
of the integrals in ()5.59p as C — ^ can be calculated, e.g. using Mathematica, and is given by 

l\ = log C + Y + 0(C), as c ^ 0, c G c± 



■dt 



logC-f + ^(C)> as C ^ 0,C G /U//U/y, 



it-C^'^ llogC + ¥ + 0(C) asC^O,CG///. 



Here I, II, III, IV denote the open quadrants of the complex (^-plane. Plugging in these asymp- 
totics in ()5.59p we get 



' 1 

log/(C)= < 



2(2^. - 1) logC + f (1 - 2v) + h{0 + 0{C), as c ^ 0, C e / u //, 



i(5 - 2v) log C + g(7 + 2v) + h{C) + 0(C), as C ^ 0, C G 
U(5 - 2u) log C + f (1 - 2v) + h{0 + 0{C), as C ^ 0, C e IV. 

In terms of the original functions /i,/2 (assuming that h is analytic) this precisely turns into 
(|5.35|) - (|5.36p for Cases I and IV. This concludes the proof of the asymptotic behavior of the 
Szego function around the origin. 

Summarizing, we have now proved that f{q) := X{q) solves Boundary Value Problem 15. 13|, 
where in (a) we only established the upper bound \ f{q)\ < C2 < +00 on Tl \ B. To obtain 
the lower bound < Ci < \f{q)\ on TZ\B, let us first define the function f{q) on TZ which is 
the solution to the same boundary value problem as f{q) except that each of the multiplicative 
jumps G{t) is replaced by its inverse 1/G{t) and the asymptotic behavior at is inverted as well. 
Proceeding as above we find such a solution f{q) which is bounded, i.e., \f{q)\ < C2 < +00 on 
1Z\B. But then the product function f{q)f{q) is analytic and bounded on TZ so by Liouville's 
theorem it must be a constant. Clearly it cannot be identically zero so it is a nonzero constant 
C 7^ 0. This shows that f{q) = C/f{q) and from the boundedness of / we then obtain the 
desired lower bound < Ci < |/(g)| with Ci = C/C2- □ 



5.5 Local parametrices near the nonzero branch points 

Near each of the branch points in {aj,bj \ j = 1, . . . , N} U {— C2, C3} \ {0}, a local parametrix 
S^"^ can be built in the standard way with the help of Airy functions, see e.g. [23]. We omit 
the details here. 



5.6 Local parametrix near the origin 

In this section we construct the local parametrix S^^^ near the origin. We will show that the 
local RH problem can be reduced to the RH problem in [T71 Section 5.6]. As a first step we 
perform a preliminary transformation on RH problem 15.101 



44 



We define the 4x4 matrix-valued function P by 

P{z) = S{-z) diag (^e±^"*/^ e^^"*/2^ e±^"*/2^ e^^"*/2^ J, for ±Im z > 0, (5.60) 

where 



J 



/O 1\ 

10 

10 

VI 0/ 



(5.61) 



Let us denote by Li = —Li, i = 1,2,3, the lens around \jf=i{—bj,—aj), (c2,+oo) and 
(— oo,— C3), respectively. We also introduce 

Xi{z) = X5-ii-z), ^ = 1, 2, 3, 4. (5.62) 

Then P solves the following RH problem. 

RH problem 5.18. The matrix-valued function P satisfies the following conditions. 

(1) P{z) is analytic for z G C\T,p, where Hp is the contour consisting of the real axis and 
the lips of the lenses L^, i = 1,2, 3. 

(2) For z G Sp, P has a jump 

P+{z)=P^{z) 

where 



diag (1, (Jp)2 , 1) , for z in M"*" and the lips of L2, 
diag {{Jp)i , iJp)^) ) for ^ in and the lips of Li, L3, 



{Jp)i = { 



_gn(Ai,_-A2,+) gi^TT 

i\ 
-1 oj • 

1 

1 _ gibz/7rign(A2 - Ai ) 

1 



.J , on (-C3,0), 

on (-00, -C3), 
on the lips of L\, 
, on the upper /lower lip of L3, 



{JP)2 = < 



I g«(A3,--A2, + ) 
1 
1^ 



-1 OJ ' 

gn(A2-A3) J I ' 



and 



( ' ^ 1 

-1 



^—V'Ki—2'Kina.j 







_g"(A3,- — A4_+) 
\ _g±!^7rign(A4-A3) \ 

1 ) ' 



on (0,C2), 

on (c2,oo), 

on the lips of L2, 

on {-bj,-aj), j = 1,...,N, 
on {-aj+i,-bj), j = 0,...,N, 

on the upper/lower lip of L\, 
on the lips of L3. 
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(3) If z ^ outside the lenses that end in 0, we have 



P{z) diagdzl'^/^, Izj-'^/^, |z|'^/2, |z|-'^/2) = ifu>{), 
P(z)diag((log|z|)-i,l,(log|z|)-i,l) =0(1), ifu = Q, 
P{z) = 0{\z\''l'^), P-i(z) = Odzl'^/^), i/-l<i/<0. 



(5.63) 



We can apply the transformation ()5.60p to the gfobal parametrix S"*^""-* as weh. That is, we 
define S^"^) ^ by 



p(°^)(^) = 5(°^)(-z)diag l^e^'^^*/^^ e^'^^^/^^ e^^^*/^ e^'^^*/^^ 



J, 



for iblmz > 0, 



(5.64) 



with J as in (|5.60p . Then the jumps for equal the jumps for P above, but with all the 

exponentially decaying entries of the form e"^'^^"^*^ removed. The behavior of near the 

origin follows trivially from ()5.47p . 

Now we observe that the jumps in the RH problems for P and are reminiscent of 

those in [171 Section 5.6], with the variable a in the latter paper playing the role of our v. Our 
construction of the local parametrix P^'^^ will be inspired by |17j . 

In the construction we have to make a case distinction between < and v > ^. To that 
end we define l;/<o = 1 if < and if > 0. The local parametrix P^'^'> is defined in a fixed 
but sufficiently small disk Z?(0, 5) around the origin, with radius 5 > 0. It satisfies the following 
RH problem. 



RH problem 5.19. We look for a 4 x 4 matrix-valued function P^^\z) : D(0, (^) \ Yip 
satisfying the following conditions. 

(1) P(°)(z) is analytic for z G -0(0,5) \ Sp. 

(2) For z G D(0, (5) n Sp, p(o) has a jump 



■^4x4 



P^^\z){ 



diag ((Jp{o));^ , ( Jp(o))3) , for z in (—(^,0) and the lips 

of Li, L3, 

diag (1, ( Jp(o))2 1 1) , for z in (0, 6) and the lips of 

L2, 



where 



p{o))l 



y<oe 



1 
-1 



_g±i/7rign(A2-Ai) 
1 



on (-(5,0) in Cases III,IV,V, 

on (—6,0) in Cases I, II, 

on D{0,5)ri the lips of Li, 

on D{0,6)ri the upper/lower lip of L3, 



(<^p(0))2 



1 1, 



ra(A3,_-A2,+) 



1 

r 
-1 0, 

1 _ 0' 

gn(A2-A3) ^ 



, on (0,6) in Cases I, II, IV, V, 
on (0, 6) in Case III, 
on D{0,5)ri the lips of L2, 
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and 



p(0) j3 



1 
-1 



-li/<oe 



1 





o±i/7rjgn(A4-A3) 
1 



on (—5,0) m Cases I,IV, 
on (—5,0) in Cases II, III, V, 

on D(0, 5)n the upper/lower lip of Li, 
on D(0,(5)n the lips of L^. 



(3) // z — ?• outside the lenses that end in 0, we have 

' p(0)(z)diag(|z|'^/M2|-''/Mz|'^/Mz|-'^/2) ^ 0(1)^ ifv>^, 

< p(0)(z)diag((log|z|)-i,l,(log|zrM) = O(l), »/z. = 0, 

p(o)(^) = 0(|z|^/2), (P(0))"^(z) = Od^l''/^), 2/-l<l/<0. 



(5.65) 



(4) On i/ie boundary of D{0,6) we have the uniform estimate 

p(o)(z) = pM(z)(/ + 0(l/n)), n — )■ oo. 



(5.66) 



Note that the exponentiahy small entries in the jump matrices on (—5,0) U (0,5) in the 
above RH problem are only present if < 0; we neglect them if > 0. The reason for this case 
distinction between i' < and z/ > is explained in [17] (with there a playing the role of our 
ly); see also the estimates in Section [5.71 below. 

To construct P^^\ we need the model RH problem for the modified Bessel function. 

RH problem 5.20. Denoting with 'jj, j = 1,2,3 the complex rays G C \ arg(^ = (j + l)7r/3}, 
we look for a 2 x 2 matrix-valued function ^^'^^^^^ such that 

(1) ^Bessel analytic inC\ U5=i 7j. 

(2) With the rays jj, j = 1,2,3 all oriented towards the origin, xj/^esscl ^^^^ jumps 



^Bessel ^Bessel ^ , 



1 0\ 






on 7i 






-1 i)' 


on 72 


1 0\ 




e-"""' 1 / ' 


on 73 



(3) Uniformly for ( ^ oo we have 



1 A i 



-1/2^ 2C1/2 



0-3 



(5.67) 



with (T3 = diag(l, —1). 
(4) As C ^ in \ argd < 2-k/3 we have 



Bessel 



(C) 



1 log Id 

1 log Id 



o 



ifu>0, 

tfu = 0, 
ifuKO. 



(5.68) 
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This RH problem has u unique solution which is given in terms of modified Bessel and Hankel 
functions, see [53]. 

We are now ready to construct P^^^ case by case, following the lines in [TTj. First we give 
the construction for Case I. We may assume without loss of generality that the lips of Li and 
L3 coincide within D{0,6). If 1/ > 0, we consider the functions 



(/.i(z) := (^Ai(z) - A2(z)±— J , 03(^) := (A3(2) - A4(^) ± 27rij , ±lmz>0. (5.69) 

It follows from (|5.62p and Lemma [5. 61 that these functions have analytic continuations to D{0, 6) 
that give conformal maps from a neighborhood of the origin onto itself, such that (l)i{x), i = 1,3, 
is real and positive for x £ (0,5). We deform the lips of Li (L3) near such that (pi maps the 
upper and lower lips of Li (L3) to the rays with angles 27r/3 and —2tt/3, respectively. 
We now define 

P(0)(.) = E(.)d.ag (a,*Bessel ^^^^ ,„,,^Bessel ^I^^ 

X diag(cJ3et(^i(")-^2W)<x3^^3gf(A3{.)-A4(.))<x3)^ (5_70) 

where cJi = ^ ^ and = diag(l, —1), and where the prefactor E{z) is analytic in D{0, 5) and 

is chosen to satisfy the matching condition on dD(0,6), see below. We use the hat superscript 
to emphasize that we are in the situation z/ > 0. Note that P^'^^ essentially decouples into two 
blocks containing the model RH problem for the modified Bessel function. With this definition, 
and assuming n = mod 3, the items (1), (2), and (3) in the RH problem for P^^^ are satisfied, 
by virtue of items (1), (2), and (4) in the RH problem for vj/Bessei^ 

To achieve the matching condition in item (4) of the RH problem for P^^\z), we take E{z) 
in ([5:7(1 as 

i?(z) = p(~)(z)(-l)-diag(a3,a3)diag(i=(_^^ "^^ , _L "^^^ 

X diag (^(^^cpy\z)y'''^\[^i,/\z)y''''^'^ . (5.71) 

Obviously E(z) is analytic for z G D{0,5) \ . Moreover, one checks that E{z) is also analytic 
across (—5,0). Finally, since (f>i,3{z)^/'^ = 0{z^^'^) as z — )• it then follows from ()5.47p and ()5.64p 
that 

E{z) = 0{z-^^^), z^O, (5.72) 

so E{z) cannot have a pole at zero. We conclude that E{z) is analytic in the disk D(0,S). By 
virtue of (|5.67p . the matching condition (j5.66p in the RH problem for P^^^ is satisfied. 

If — 1 < < 0, we cannot simply ignore the jumps on the real axis. The local parametrix 
P^^^ is then constructed in the following form 

P(0)(z) = pW(z)diag(e-"^i(^),e-"^2(^)^g-nA3(^)^g-nA4(^))g(^) 

X diag(e"^i(^\e"^2(^)^gnA3(^)^gnA4(^))^ (5_73) 

where P^*^^ is the parametrix for the case z/ > given in (j5.70p . and Q{z) is a piecewise constant 
matrix. More precisely, following the idea in |17[ Section 5.6.3], we have 

Q(z) = I - —-—-E2,3, (5.74) 
2zsm(7rzvj 
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for z in the region bounded by (0,5) and the upper Hp of Li, 



Q{z) = l 



^2,3, 



2i sin(7ri^) 

for z in the region bounded by (0, 5) and the lower lip of Li 



(5.75) 



l\ 1 e""' \ 

2i sin(7ri/) 2i sin(7rf) ' 

1 


Vo 



2i sin(7r!^) 
1 





2i sin{7n^) 



1 



(5.76) 



for z in the region bounded by (—5, 0) and the upper lip of L 



/I 

1 



Vo 



2jsin(7n/) 2isin(7rj/) 

2jsin(7r;/) 2isin{7r!^) 

1 

1 / 



(5.77) 



for z in the region bounded by (—5,0) and the lower lip of L\. With Q given in (|5.74p ~ ()5.77p . 
one can check that P^^'^ defined in (|5.73p indeed satisfies the items (l)-(4) in RH problem l5.19l if 
— 1 < < 0. For the jump condition (2) we use that n = mod 3. For the matching condition 
(4) we also need the inequalities Re A3^4(z) < Re Ai^2(-2) for z in a neighborhood of the origin, see 
Lemma [5. 41 and (|5.12p ~ ()5.13p . Moreover, this construction actually works as long as sin(7rz^) ^ 0, 
i.e., V ^ NU {0}. 

For Case II, the RH problem for P^^^ is exactly the same as the one considered in ]yi\ 
Section 5.6]. Also the construction of P'^^^ in Case IV is similar to Case II. We thus omit the 
details for these cases. 

For Case III, the Bessel parametrix will appear in the middle block. More precisely, by 
setting 



4>^{z) := [ Xsiz) - \2{z) ± ^) = -clz + 0(z2), 



± Imz > 0, 



(5.78) 



we have 



PW(Z) = E(z)diag (^^--/2,vI;Bessel {^^^^^^ .Z^l^ 

X diag (l,a3e±"f^'^3et(A.(.)-Aa(.)).3^ 1^ 
and the analytic prefactor E is given by 



±lmz>0, (5.79) 



E{z) = P(°°)(z)diag (l,(j^e^'^''\l 



X diag ( z' 



^2 ) ' ^ 



±lmz>0. (5.80) 



This describes the parametrix if > 0. If — 1 < < 0, we define P^^\z) by ([5:73]) . ([5779]) . where 
now Q{z) is a piecewise constant matrix given by 



Q{z) = I + 



2 sm.{'Kv) 



{^2,1 + -£-4,3) 



(5.81) 
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for z in the region outside the lens, 

/ 1 



1 sin(7r!^) 

i 

2 sin(7ri^) 





V 



0\ 

1 
10 

i i 1 1 

2 sin(7rj^) 2 sin(7n/) / 



(5.82) 



for z in the region bounded by (0, 6) and the upper hp of L2, and 

/ 1 

Qiz) = 



2 sin{7ri^) 



V 



2 sin(nu) 







1 





2sin(7rj^) 2sin(7n/) / 



o\ 




1 



(5.83) 



for z in the region bounded by (0,(5) and the lower hp of L2. To check the matching condition 
(4) in RH problem 15.191 we need the inequalities ReAi(z) < ReX2,3{z) < ReX^^z) for z in a 
neighborhood of the origin, see Lemma 15.41 and (|5.14p . 

Finally we build the local parametrix in Case V. The construction is much simpler in this 
case. We set 

p(o)(z) =p(°^)(z). (5.84) 

This describes the parametrix if u > 0. If — 1 < < 0, we define p(°)(z) by (|5T3]l . (fTOD . where 
now 

/ 1 0\ 

1 



Q{z) 



2 sin{7n^) 






2 sin(7rj/) 

1 

1 



(5.85) 



2 sin{7ri^) 

if iblmz > 0. One easily checks that conditions (l)-(4) in RH problem 15.191 are satisfied. 

Finally, by tracing back the transformation (j5.60p . we obtain the local approximation S^^^ 
for 5 from 



5.7 Final transformation 

Using the local parametrices described in Sections 15.51 and [5.61 and the global parametrix S^°°\ 
we define the final transformation as follows 

{S{z){S^^'^^ (z))^^ , in the disks around {aj,bj,ci, —C2,cs} \ {0}, 
5(z)(5(°)(z))-\ in the disk P)(0, 5) around the origin, (5.86) 
S{z){S^'^\z))-\ elsewhere, 

where S, S^'^y, S'^^\ and S'^^) are respectively defined in (ICTj) : Section I^M), (I^TTOj) . 
(iLZal), (I57r9l> . and (lOTD : and (IQHjl . 

From our construction of the parametrices, it follows that R satisfies the following RH 
problem. 

(1) R is analytic in C \ S^j, where the contour S/j depends on whether < or > and is 
different for all five cases. 

(2) R has jumps R+ = R-Jr on that satisfy 

jR{z)=I + 0{l/n), (5.87) 
uniformly for z on the boundaries of the disks; 

Jji(x) =1 + 0{x^e-^"), (5.88) 
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on (-(5,0) (Case I), or on (0,6) (Case III), or on {-5,0) U (0,5) (Cases II, IV and V), for 
some constant c > 0, if > 0; and 

J,j(z) = / + 0(e-^"l^l), (5.89) 

on the other parts of S/j, for some constant c > 0. 

(3) R{z) =1 + 0{l/z) as z^oo. 

In case > 0, the estimate ()5.88p is not trivial. We prove it in Case I. In that case R{z) is 
clearly analytic in D(0, 5) \ (—(5,0] with the following jump on (—(5,0) 

R-{x)-^R+{x) = p|°^(-x)Jp(-x)-ipl°^(-x)-\ 

= Pf\-x) (l - e«(A3,-A2,H-){-x)^2 3^ p'^\-x)-\ 

= / _ e"(^3.— ^2.+)(-^')p(0)(_2;)^2,3P^°H-a;)~^ 

Now for z/ > the matrix P^^\-x)E2,zP^^\-x)-^ = 0{\ xY) as X — )• 0, X < 0. Indeed (following 
the proof of [1% Lemma 5.6]) we observe that 

P^'^\-x)E2,zP^'^\-x)-^ = P^^\-x) {0 1 0)^(0 1 0)p(°)(-x)-i. 
Then if z^ > we find from (|5.65|) and the fact that det P^^^ {z) = 1 that both factors 

p(°)(-x)(0 1 0)^ and (O 1 O) P^^\-x)-^ 
behave like 0(|x|^/2) 

as X — >• 0, X < 0, which proves the statement for > 0. In case = 0, this 
approach does not work as it would lead to a bound ©(log |x|). However, we do find that the 
first factor pW(-x) (0 1 O) remains bounded as x — )• 0, x < 0. For the second factor we 
observe from ([5^ and det ^^^'''''''(C) = 1 that 

(^B„,)-'(,,_(0(log|C|) 0(k,g|C|)^_ ^^^^ 

Combined with (j5.70p and the fact that E{z) is bounded near the origin we obtain that 
(O 1 O) P*^''''(— x)~^ is also bounded as x — >• 0, x < 0, which proves the statement in 
case V = 0. Cases II, III, IV, and V can be similarly treated. 
Then from standard arguments we may conclude that 

as n — )• oo, uniformly for z in the complex plane outside of S^. 
5.8 Proof of Theorem [276] 

We follow the approach in [231 Section 9.3] and start the proof with a lemma. 
Lemma 5.21. For every x ^ S{^i)\\Ji^^-^{D {at-, e) ^ D{hk,e)) we have 

S-^\y)S^{x)={^^ + ''[^-y^ *), asy^x, 

uniformly in n. The * entries denote unimportant 2x2 blocks. 

Proof. The proof is standard, see e.g. [231 Lemma 9.9]. D 
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Recall the formula ()2.1ip that expresses Kn in terms of Y. The idea is then to write this 
expression in terms of R instead of Y by applying all respective transformations 

introduced in the steepest descent analysis. Meanwhile, we let n tend to infinity, so that we can 
exploit the conclusion of the steepest descent analysis ()5.90p . 

Let x,y > 0. First, unfolding the transformation Y X given in (|5.4p we get 

Kn{x,y) = —— (0 Wo,n{y) Wi^niy) W2,n{y)) 

2tti[x — y) 

xdiag(l,Z)(z)-iTy-^(y)e"Q+(^))x-i(y)X+(x)(l O)^. 
Using ()5.3p this boils down to 



u/2 -u/2 n{ei^+{,y)~V{y)) 

Kn{x,y) = ^- ^-Kiix-y) ^ ° X-\y)X+{x) {l O)^. (5.91) 

Next, by the transformation X ^ U described by ()5.19p we obtain 

„n(A2,+ (y) — Ai,+ (a;)) 

^n(y(y)-v{x))j^^^^^y^^yU/2^-u/2^^^ ^ q U^\y)U+{x) [l 0)^. 

2tti[x — y) 

The opening of global lenses [/ i— )■ T in (|5.28p does not effect the above expression. The 
opening of the local lens T i— )• 5 in (|5.29p . however, does have impact 

„,iy/2^-iy/2 

e<ny)-n^))Kr,{x,y) y 



2'Ki{x — y) 

where Xs(p.i) denotes the characteristic function of the set S{fii). 

Let X G S{fii) \ {ak,bk \k = l,...,N}. The factors e"(^(?')-^(^)) and y^'/^x-"/^ disappear 
as y — >• X. Then Lemma 15.211 vields 

gri(A2,+ (y)-A2,+ (a;)) _ gn(Ai,+ (?/)-Ai,+ (a;)) 



Hence, 



1 Id 
lim -Kn{x,x) = — — (Ai,+ (x) - A2,+(2;)) 
n-s>oo n Ztti ax 



'X 

where the second equality follows from (j5.15p . the third is taken from |231 p. 116], and the last 
one is ^TM . 
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6 Triple scaling limit in the quadrat ic/linear case 



In this part we study the very concrete case of the chiral two-matrix model with potentials 

V{x)=x, W{y) = ^+ay. 

For this case we were able to construct a phase diagram in the (a, r)-plane, see Figure [TJ Very 
remarkable is the occurrence of a multi-critical point for the parameter values a = — l,r = l. In 
this part we will study a triple scaling limit to this point leading to the chiral version of the main 
result in [20]. An essential point in the proof is the construction of the local parametrix around 
zero. In this construction we will make use of the solution to a model RH problem introduced 
in dS]. 

The proof of Theorem 12.81 is again based on a a steepest descent analysis. This analysis is 
very similar to the one performed for the non-critical cases. In fact the first four transformations 
are almost exactly equal. The only difference is that we will need to modify the A-functions, 
which will be introduced first. 

6.1 Modified A-functions 

Let us first introduce an auxiliary parameter 7 that is completely determined by a and r but 
will prove to be convenient for notation. We define 7 = 7(0, r) as the solution of 

0x2/3 = ^-972 + 5x^/37, (6.1) 
7 

that tends to 1 as r — )• 1 and a — )• — 1. In the triple scaling limit, i.e. we let a and r depend on 
n as in (|2.25p while n — )• 00, we have 

7 = 1 + ian-V3 + ^u_y ^ 47^) „-2/3 + (g,2) 

Lemma 6.1. There exist functions Xj, j = 1, 2, 3, 4, analytic on C \ M that satisfy the following 
conditions. 

(a) As z ^ 00 we have 

Xi{z) = z-log{z) + ii + 0{z-^) , (6.3) 
\2{z) = 6i{z) + ^ log z + £2 + Cz-i/3 + Dz-2/3 + o (z-i) , (6.4) 
^ , , r 1 log z + ^3 + Ca;z-V3 + 1)^2^-2/3 + Q (^-1) c+, , , 

A3 (z) = 6*2(2;) + < ? „ -, ; -,; (6.5) 

^' ^' \\\ogz + U + Cu^z'^/^ + Dujz-^/^ + 0{z-^) inC-, 

(^logz + U + Cu'z-^/^ + Dujz-y^ + 0{z-^) znC+, 

A4(2;) = 6^3(2;) + < ■? , ,„ „ „,„ ; (6.6) 

I i log z + £3 + ^6^2-1/3 + ^^2^-2/3 ^0(^-1^ ^ ^ 

where C,D gM. and ij, j = 1,2,3,4, are constants satisfying 

l^-l2 = l2-U = \^i. (6.7) 

(b) There exists a constant c > such that the modified X-functions satisfy the following jump 
conditions. 

(i) On R"*" we have 

(6.8) 
(6.9) 
(6.10) 





— A2,=p 


on 


(0,c), 




= 


on 


(c,oo), j = l,2. 


^3,± 


= A4,=p 


on 


(0,00). 
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(ii) On M we have 



\j,+ = Xj,--27ri, j = l,4, (6.11) 

A2,± = As,^ ± 27ri. (6.12) 

(c) In a neighborhood of the origin we have 

F{z)z'^l^ + G(z)zi/2 + H{z)z^/^ in C+, 

iF{z)z'^/^ - G{z)z^l^ - iH{z)z^/^ + 2m in C_, 

'iF{z)z'^/^ - G{z)z^/^ - iH{z)z^/^ + 2m in C+, 

F{z)z'^/* + G{z)z^l'^ + H{z)z^l^ inC, 

-iF{z)z^l^ - G{z)z'^/^ + iH{z)z^/^ + 2m in C+, 

-F{z)z^/^ + G(z)zi/2 - H{z)z^'^ in C_ , 

-F{z)z^l^ + G{z)z^''^ - H{z)z^l^ m C+, 

-iF{z)z^l^ - G{z)z'^/^ + iH{z)z^/^ + 2m in C_ . 

where F, G, H, K are analytic functions satisfying 

F{0) = 4e3-/S'/' (-27' + ^ + r^/^) , (6-13) 

G(0) = 2.7-^/^ (§7' - ^ - fr^/^) , (6.14) 

//(O) = 2e-/S-^/^ (17^ - ifc + ir^/^7) • (6.15) 




Proof. Define 



Xj{z) = 2X^^{^), i = 1,2,3,4, z£C\R, (6.16) 



where A^^ denote the modified A-functions introduced in [2Q^ Section 3.4]. Then (a) follows from 
|20[ Lemma 3.11]. (b)(i) is direct from |20[ Lemma 3.10]. To prove (ii) we need some symmetry 
conditions on A^c. By [201 (3.15), (3.23)^(3.24), (3.29)-(3.30)], it follows that for x > 

Xf^iix) = X^^{-ix) - TTi, j = 1, 4, 

Xf^^iix) = X^^i-ix) + TTi, j = 2, 3. 

This, together with ()6.16p . implies (b)(ii), where we take c = (c'^^)^. Finally (c) follows from 
\2Q\ Lemma 3.12] where we put 



>(z) = 2FNC(^), Giz) = 2GNC(^), 
Hiz) = 2FNC(^), k{z) = 2i^NC(^)^ 



z G c\: 



□ 



Remark 6.2. The constants 7, c, etc. and the functions F, G, H, K, Xj, etc. all implicitly depend 
on a and r (and thus also on n via the triple scaling limit). When dealing with these functions 
associated with the critical values of the parameters a = — l,r = l, we add a star to the 
notation. Thus, we write -f* , c* , F* ,G* , H* , K* , X* . . . 

6.2 The transformations Y X U i-j-Th-S' 

The first transformations yi— t-Xi— T-t/i— )-Ti— )-5of the steepest descent analysis in the critical 
case are almost the same as before. The main difference is that we use the modified lambda 
functions rather than the ordinary ones. Apart from that also some details have to be changed. 
We list them here. 
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Transformation Y X. 

This transformation is exactly as in Section 15. 1[ 



Transformation X U. 

Definition 15. 71 of U{z) has to be slightly changed. Besides the fact that we replace the A-functions 
by the modified A-functions, we also have to change the definition of the diagonal matrix L by 

^ f-diag(£i,£2,4,4) forlmz>0, 

I -diag(£i, £2,4,^3) forlmz<0. 

Then, under the assumption that n = mod 3, U solves RH problem 15.81 if we put := 1, 
ai := 0, bi := c, and C2 = C3 := 0. 

Transformations U T S: opening of lenses. 

Here, we open the local lens Li around [0, c] and unbounded lenses L2 around and L3 around 
M"*". We want to do this such that the off-diagonal entries of the jump matrices on the lips of 
these lenses tend exponentially fast to zero as n — )• 00. This is only possible outside a shrinking 
disk around the origin. It will be sufficient for our purposes to let this disk shrink with speed 

D (o,pn-^^A = {z G C I l^l < pn-^/^}, 



where /) > is a constant that will be chosen sufficiently small later on. 
Lemma 6.3. The lenses Lj, j = 1,2,3, can be opened (independently of n) such that 

Re(Ai(z) — A2(^;)) < —dn^^^"^, for z on the lips of L\ but (6.18) 

outside D (o,pn^^^^ 



Re(A3(z) - A2(z)) < -(in~^/2max(l,|z|2/3), for z on the lips of L2 but (6.19) 

outside D (o,pn^^^^ 



Re(A3(z) - A4(z)) < -dn"^/^max(l,|z|2/^), for z on the lips of but (6.20) 

outside D (o,pn^^^^ 



for sufficiently large n and a fixed constant d > 0. Moreover, there exists a constant d' such that 

Re(A2(x) - Ai(x)) < -d'x, x E (c* e, 00), (6.21) 
where e > is a small number. 

Proof. Estimates (|6.18p ~ ()6.20p are immediate from ()6.16p and \20\ Lemmas 4.6 and 4.9] where 
we define Lj = (^LJ'^^ . Estimate ()6.2ip follows in the same way from the following result that 

holds in the context of |20j but was not explicitly mentioned there: there exists d > such that 
for sufficiently large n 

Re (A^C(x) - Af C(x)) < -dx\ x e (-00, -(c*)^^ - e) U ((c*)^^ + e, 00) . 

□ 

In this way we arrive at the following RH problem for 5. 

RH problem 6.4. The matrix-valued function S is the unique solution of the following RH 
problem 
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(1) S{z) is analytic for z G C \ E5, where is the contour consisting of the real axis and the 
lips of the lenses Li, i = 1,2,3. These lenses are chosen such that the estimates in Lemma 
fOI hold. 



(2) For z G Ss, S has a jump 

o ( \ = Q ( )i "^^^S {{Js)i (z), iJs)^ (z)) , for z in M+ and the lips of Li, L3, 
^+{Z) ^-{.z) <^ ^.^^ j^^^i^ ^j^^^ (z),e-'^^^) , for z in R~ and the lips of L2, 



where 



iJs)i = { 



I g"(A2,+-Al,-) ^ 

,0 1 

1 0\ ' 

gn(Ai-A2) 2 / ' 



iJs)2 = { 



and 



1 
-1 Oyl' 

I _g±i''rign(A3-A2)^ 
1 



.-1 »)' 

2, 

' 1 0^ 

gn{A3— A4) ^ 



{Js)3 = < 



m 

on (0, c), 

j , on (c,oo), 

on i/ie lips of Li, 
on the lips of L3, 

on (—00, 0), 

on the upper/lower lip of L2, 

on (0, +00), 
on the lips of Li, 
, on the lips of L3. 



(3) As 2; — )• 00 with iblmz > 0, we have 

S{z) =[l + (z-i)] diag (1, 2^/3^ 1) diag (z'^/^^ 2"^/^^^) diag (l, 1, a^,a^^) . 



(4) S{z) has the same behavior near the origin as X{z), see ()5.6p . provided that z — )• outside 
the lenses that end in 0. 

We will construct local and global parametrices for P{z) as described in Section [5161 i.e. P{z) 
is established from the solution S{z) to RH problem 16.41 bv the transformation ()5.60p . More 
precisely, we have 

RH problem 6.5. The matrix-valued function P satisfies the following RH problem. 

(1) P{z) is analytic for z G C \ Sp, where T,p is the contour consisting of the real axis and 
the lips of the lenses Li = —Li, i = 1,2, 3. 

(2) For z £ Tip, P has a jump 

diag (1, ( Jp(z))2 , 1) , for z in and the lips of L2, 



P+{z) = P^{z) 



diag{{Jp{z))^ ,{ Jp{z))^) , for z inM and the lips of Li, L3, 
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where 



PI 



^ ^ 1\ 

-1 Oy) ' 

I \ _g±!^7rign(A2-Ai) 

\0 1 



{JP)2 = { 



1 
-1 Oyl' 
1 0^ 

gn(A2-A3) 



and 



_g'^{A3,-— A4,+) 

\ _g±i^7ri+n.(A4-A3)\ 

.0 1 j 



U2, 



on {—CO, 0), 

on the lips of Li, 

on the upper /lower lip of L^, 

on (0, 00), 
, on the lips of L2, 

on (— c, 0), 
on (—00, —c), 

on the upper /lower lip of Li, 
on the lips of L3, 



where 



5.22) 



(3) If z ^ outside the lenses that end in 0, we have 

" P{z)dmgi\z\''/^,\z\-''/'^,\z\''/^,\z\-''/'^) = 0(1), ifu>0, 
P{z)dmg{{log\z\)-\l,{log\z\)'\l)=0{l), ifu = 0, 
P{z) = Oilzl"/^), P-'^{z) = 0{\z\''/^), if-l<iy<0. 

6.3 Global parametrix 



(6.23) 



By ignoring all the exponentially small terms in the RH problem for P, we obtain a RH problem 
for p(°°) . This global parametrix can be constructed by first constructing S^°°^ along the 

lines of Section [5.41 and then again applying the transformation 5"^°°) 1— )• as given by (j5.64p . 
Its behavior around the origin is given by 



as z — >• 0, zb Im z > 0, 



for constant matrices P± , Q± . Clearly 

P±Q± = Q±P± = 0. 

Furthermore, in view of the jump conditions for it is readily seen that 



P+ = P_ diag ( 1, ( J 1 , 1 



P+ = P-e^'^^/^diag 



1 
-1 



1 
-1 



(6.24) 

(6.25) 

(6.26) 
(6.27) 
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Eliminating P+ from these two formulas we get 



Iterating this relation we see 












-1 


°\ 




-1 
























1 






Vo 


-1 









/ 









= P_ 







—i 









i 












V 


-i 





0^ 





(6.28) 



Similarly, we have 



Q+ = diag(l,( J Q^),l)g- 



-n /o -r 



/O -i\ 
i 
-i 



(6.29) 
(6.30) 

(6.31) 



\i / 
These relations will be helpful later. 

6.4 Local parametrix at c 

As in Section 15.51 the local parametrix S^^^ around c can be built in the standard way with the 
help of Airy functions, see e.g. [23]. We omit the details here. 

6.5 Local parametrix at the origin 

In this section we build a local parametrix P^^^ near the origin. Here the analysis is essentially 
different from the noncritical situation discussed in the previous section. 

Transformation of the RH problem for M(C) 

To build the local parametrix near the origin we will use a slightly modified version of RH 
problem [2.71 We put u = u + 1/2 and set 



iV(C) = diag( (J M,l,zjM(C)diag((^. 



(6.32) 



The jumps for N are shown in Figure [TUl 

The asymptotics of as (" — ^ c« is given by 



mc) = ( ^ + ^ + ^ + o (^) ) diag(c-^/^ (-c)-^/^ (-c)^/^ e^') 



Ni N2 

c 



X ^diag (e-^i(^)-*^,e-^2(f)+*^,e'^2(«+*^,e'^i(^)-*^ ) (6.33) 



with 



A 



V2 



/I i\ 

1 -i 

-i 1 

\i 1/ 



(6.34) 
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Figure 10: The figure shows the jump matrices J^, A; = 0, . . . , 9, in the RH problem for A'^ = N{C,). 



and the behavior of A^(C) for C — t- is given by 

iV(C) = 0(n, N-\0 = O{C), ifi><0, 

and 

iV(C) diag(C, C', r , C-') = 0(1), C G f^i U 1^8, 
A^(C)diag(C-^r,C-^n = 0(1), CGf^sUl^e, ifP>0. 

iV(C) = 0(c-^), C0(f^iUi^3Ui^6Ui^8), 

Note that the behavior at infinity can be rewritten as 

N{c) = diag(c-^/^ (-c)-^/^ (-c)^/^ e'')A + ^ + ^ + ^ (c^) ) 



X diag (e-^i(«-*^, e-'^2(f)+*^, e^^CO+tf^ e^i(C)-tC j (6.35) 
as C — >• 00 within the upper/lower half plane. Here 

iVi,± = diag {l, e^^'l^, 0, o) A^i diag (^0, 0, e^'^*/^ l) I. (6.36) 
For further use, we record the symmetry relation (see also [16j) 



iV(-C;ri,r2,s,t) = diag 



1\ ( 1 
-1 O/'l-l 



xiV(C;r2,ri,5,i)diag( (J ^]]. (6.37) 



Note that the order of ri and r2 differs at both sides of the equality. 
A corollary of this relation and (|6.35|) is 

^i,T(n,r2,s,t) //O 1\ / 1 

— diag ' ' ' ' 
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Construction of the local parametrix 

Here we construct the local parametrix P^^^ around the origin. We will do this inside a shrinking 
disk D (0,/9n~^/^), i.e. we want to solve the following RH problem. 

RH problem 6.6. We look for P^'^^ satisfying the following conditions. 

(1) P^'^{z) is analytic for z £ D{0, pn^^^^)\I]p, where D{0, pn^^^^) denotes the disk of radius 
pji-^/'i around 0. 

(2) satisfies the jumps 

pjO) ^ j^^ p ^(Q^ ^-1/3)^ 

where Jp is the jump matrix in RH problem \6.5[ 

(3) As n — )■ oo, we have that 

p(0\z) = + Z{z) + ©(n-^/^^)) P'-°°\z), uniformly for \z\ = pn'^/^, (6.39) 
where 

Z{z) = as n ^ oo, uniformly for \z\ = prT^^'^ . (6.40) 

We define the local parametrix as follows (compare with |16j): 

P(0)(z) = En{z)N (n^/''f{z);r,{z),r2{z),n^/h{z),n'/h{z)) A{z), (6.41) 

where 

A(z) = diag (e"^i(^),...,e"^^(")) , 

and, as before, \j{z) = A5_j(— z), j = 1, 2, 3, 4. The parameters n, r2, s, and t depend in a mild 
way on z and also on n, although this is not indicated in the notation. The prefactor En{z) 
does depend on n and is analytic in a neighborhood of the origin except on the negative real 
line where it has a jump. 
We set 

f{z) = z^l\ (6.42) 

where we put the branch cut along (— oo,0]. Furthermore we define for z € L'(0, pn^^^^) \ 
(_p„-i/3^0] 

ri(z) = ^e--/^i/(-z) - y^^l\F{-z) - F{0))z-'/\ (6.43) 

r2(z) = ^e-^'/'^H{-z) + ^e'^^/^(F(-z) - F{0))z~^/\ (6.44) 

s{z) = le-3-/4p(o), (6.45) 

t{z) = iG{-z), (6.46) 

where F, G, and H are the analytic functions introduced in item (c) of Lemma 16.11 Clearly 
t{z) and s{z) are analytic in the disk D{0, pn~^/^) for n large enough. ri(z) and r2{z) are only 
analytic in D(0, pn^^/^) \ (— pn^-*^/^, 0], but they are each others analytic continuation across the 
interval {—pn^^/^,Q). 

In the construction of the parametrix, apart from the analytic structure of the parameters 
fiiz), f2{z), s{z), and t(z), we will also need their asymptotic behavior as n — ?■ oo. 
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Lemma 6.7. The functions s{z) and t{z) are analytic for z ^ D (0,pn ^/^) and n sufficiently 
large. The functions ri{z) and r2{z) are analytic in D (0,/9n~^/'^) \ {—pn~^/^,0] and satisfy 



i"i,±{^) = f2,^{x), —pn < X < 0. 
Moreover, there exists a constant G G C, independent of z and n, such that 



.47) 



lim, 



n^oo fjijT' ^ ^) — ^) 

\imn-^oD'n'^^^s{n^^/^z) = i(a^ — bb), 
lim^^oo ni/3t(n"^/3z) = 2a + Gz, 



(6.48) 



for \z\ < p and j = 1, 2. 



Proof The analytic structure of these functions for sufficiently large n was already discussed. 
The limiting behavior for rj is obvious. To obtain the limiting behavior for t we write 



'Giz) = G{0) + Giz + O{z^), 
G*{z) = G*{0) + Glz + O{z^), 



as z 



0. 



Hence 



z 



^^^GiO) -iGiz + Oin' 



1/3N 



as n — 7- oo, which follows from ()6.46p . Using this, ()6.14p . and the fact that Gi — t- as n — ?■ oo, 
we obtain 



lim n^/^t( 



n i/^z) = -2 lim ni/^-'/' - jz - jr^^^l 



iG\z. 



By inserting the limiting behavior for 7 and r as given in ()6.2p and ()2.25p . we obtain the 
statement. 

The proof for s{z) is easier. Using (j6.45p we obtain 



n^/h{n-^'''z) 



e~F(0). 



Using ()6.13p we get 



lim n^^^sin ^^'^z) 



2 lim n2/37i/4 



-27^ + h + r"'^ 



which in combination with ()6.2p and ()2.25p . finishes the proof. 



□ 



From Lemma 16.71 it follows that (|6.4ip is well-defined (postponing the definition of En for 
a moment). Indeed, it follows from standard arguments that if the solution to RH problem 12.71 
exists, it depends analytically on the parameters r^, s, and t. Combining this observation with 
Lemma 16.71 we see that for the choice of rj,s and t we made, the solution to RH problem 12.71 
exists and hence (|6.4ip is well-defined for sufficiently large n. 

Lemma 6.8. Given definitions (|6.43p - (|6.46p . the following formulas hold modulo 27ri. 

' n(\i{z) + zK{-z)) = il)i{n'^/^z^/'^;ri{z),n'^/^s{z)) +nt{z)z'^/^, 
n(X2{z) + zK{-z)) = V'2(n2/^2^/^; r2(z), n2/3s(z)) - nt{z)z'^/'^, 
niXsiz) + zK{-z)) = -V'2(n2/3zi/2. ^^(z), n^/^s{z)) - nt{z)z'^l^, ^ ' ' 

n(Xi{z) + zK{-z)) = -V'i(n2/3zi/2; ri{z), n^/^s{z)) + nt{z)z^/^. 



Proof. This is a straightforward verification. 



□ 
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It remains to define the prefactor 

X diag {n'/^f{z)'/\n'/^{-f{z))'/\n-'/'{-f{z))-'/\n-'/^f{z)-'/'') (6.50) 

with A in ()6.34p . This prefactor is not analytic in the fuh disk D (O, /9n~^/'^) , but has a jump as 
indicated in the following lemma. 

Lemma 6.9. En{z) is analytic in D{0, pn~^/^) \ with a jump 

En,+ {x) = En,^{x)dmg(^(^^^ o)'(-l o))' -pn-^/'<x<0. (6.51) 

Proof. It follows from the definition of En{z) and the fact that P^°°\z) is analytic for z £ 
—pn ^/^,pn -^/^) that En{z) is also analytic for z G D{0,pn ^/^) \ M. For 
X G (0,pn^^/^), one checks that the jumps of p(°°) and of the rightmost factor in (j6.50p cancel 
each other out so that £^(2;) is analytic for x G (0, pn~^^^). A similar calculation yields the jump 
(f63T]) of En{x) for X G {-pn'^/^, 0). □ 

Now let us check that P^^\z) in (|6.4ip has the correct jumps. This is straightforward for 
the jumps on the 5 rays which are not M~. To check the jump for x £ M~, we calculate (see also 
il6j) 

Pf (a;) = (x) diag ^ ^ ° ^ ° ^ 



1 oy'V-1 0^ 

X A^(n2/3/+(x); ri^+{x),r2,+ {x),n^/^ s{x), n^/^t{x))A+{x) 
En,-{x)N{-n^/^U{xy, r2,+ {x),ri^+{x),n^/h{x),n^/h{x)) 

-1 oi'l-i 0, 



X A_(x) diag 



.Pi»',.)aia.((_», ;),(_», I 

where in the first step we used (|6.4ip and (|6.5ip . and in the second step we used (|6.37p and 
commuted with A(x). The final equality follows from ()6.47p and ()6.4ip . This yields the required 
jumps on M~. 

The following lemma states how the local parametrix matches with the outer parametrix on 
the boundary of the shrinking disk D (0,pn~^/^), settling item (3) in the local RH problem 16.61 

Lemma 6.10. For z on the shrinking circle the following matching condition uniformly holds 
as n —7- 00 

pW(z) (p(°°)(z))"^ =/ + Z(z) + 0(n-i/i2^, \z\=pn-^/^, (6.52) 
where _ 

n'-'-^z 

and where Ni^±, P± and Q± are defined in (|6.36p and (j6.24p . respectively. 
Proof. From (j6.4ip . it follows that 

p(o)(^) ('pM(2))"^ =^„(^)iV(n2/3/(2);ri(z),r2(z),n2/3s(z),ni/3t(z))A(z) fp(~)(z))"\ 
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If \z\ = pn then \'n?/'^ z^^'^\ — )• oo as n — )• oo. Hence, we can use ()6.35p . ()6.50p and Lemma 
T8]to obtain, for itlmz > 0, 



^1 



This, together with (fOD . leads us to (|632l) and (1633]) . □ 
Note that Z(z), as defined in ()6.53p . remains bounded as n — )• oo with \z\ = pn^^f^ . 

6.6 Transformation S ^ R 

Using the global parametrix S^°°\z) and the local parametrices S^'^^ and S^^^ around c and 
we define the transformation S R as 



R 



'S{S^^^)-\ in the disk D {0,pn-^/^) around 0, 

in a fixed disk around c, (6.54) 
^S{S^^^y\ elsewhere. 



Then R is defined and analytic outside S5 and the two disks around and c, with an analytic 
continuation across those parts of on which the jumps of the parametrices coincide with those 
of S. What remains are the jumps on a contour that consists of the two circles, the part of 
the interval [c*,oo) outside the disk and the lips of local and global lenses outside the disks. 

By setting the orientations of circles to be clockwise, we have that R satisfies the following 
RH problem. 

RH problem 6.11. The matrix-valued function R satisfies the following RH problem. 

(1) R is analytic in C\ S^. 

(2) R satisfies the jump relation = R-Jj^ on T,^ with jump matrices 



'^-CO) (5.(00))-! ^/jg boundary of D {O, pn"^/''^) , 

5'(c)(^5'(oo)~)-i ^/jg boundary of the disk around c, 

^5(°°)J5(5(°°))-1 elsewhere on Sj^. 



(3) As z ^ 00, we have ^ 

R{z) = l^O{\lz). 

The jump matrix is not close to the identity matrix on the shrinking circle around 0, 
since by Lemma 16.101 we have 

J~(z) = 5(°)(z)(5(°°)(z))-i =p(°)(-z)(p(°^)(-z))-^ =/ + Z(-z) + 0(n"i/i2)^ asn^oo 

uniformly for \z\ = pn~^^^, with Z(z) = 0{1). 

The other jump matrices, however, are close to the identity matrix as n gets large. 
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6.7 Final transformation R*-^ R: nilpotent structure. 

The presence of the bounded term Z{z) in Lemma |6. 101 requires an extra transformation. 
Lemma 6.12. The function Z[z), defined in ()6.53p . has the following properties. 

(a) Z{z) is meromorphic in a neighborhood of zero with a simple pole in zero. Hence we can 
write 

where Zq = Res(Z, 0) is independent of z, and Z{z) — Zq/z is analytic in z. 

(b) Z{z) is nilpotent of degree two, moreover 

Z{zi)Z{z2) = 0, for any zi,Z2 in a neighborhood of the origin. 

Proof. To prove (a) note that it is clear from (|6.53|) that Z{z) is analytic in a neighborhood of 
zero with cut along the real line. It is then sufficient to show that there is no jump on the real 
line. When x G (0, p), there are no jumps for the functions ri, r2, s, and t. Hence, it follows 
from dOH]) and IK^ that 



1,+ 



/ 1 I 

-i -1 

I -i 

\-i 1/ 



iVi 



/I -i\ 

-i 1 

-1 -i 

\i 1 / 



(6.55) 



On account of (|6.26p and (j6.29p . we obtain from (|6.53p that 



/ 1 i\ 

1 -i 

i 1 

\-i 1/ 



iVi - 



/I 

1 i 

-i 1 

Vi 1 / 



Z^(x) 



where the second equality follows from (|6.28p and (|6.3ip . On the negative real line we also need 
(|6.38p . Indeed, for —p < x < and C± ■= n'^^^x]^'^, we have 



P+ iVi,+ (ri,+ (2;),r2,+ (x),n2/3s(x),nV3t(x) 



3/8 



1/3 1/4 

P- iVi,_(r2,+ (a;),ri,+ (x),n^/^s(x),n^/^i(x) Q 



3/8 



X 



3/8 



3/8' 



where we have used (j6.38p and the relations (|6.27p and (j6.30p . Combining this with (|6.47p and 
the fact that ^+ = — we see that Z{z) is continuous across (— p, 0) and, thus, also analytic in 
a punctured neighborhood of the origin. Recalling (|6.53p we then obtain (a), 
(b) follows from (|6.53p and the observation that Q±P± = 0, see ()6.25p . 



As a corollary we also get 

Z(z)Zo = 0, ZoZ(z) 
Then we define the transformation R as 



0, 



R{z) 



R{z) [I + Z{-z) + fi) , for z G D (O, pn" V^) \ ^-^ 
R{z) (/ + fi) , for z G C \ (Z) (0, y ^^). 



□ 



(6.56) 



(6.57) 



Then R is defined and analytic in C \ S/j where = and R satisfies a RH problem of 
the following form. 
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RH problem 6.13. 



(1) R is analytic in C \ S/j. 

(2) R satisfies the jump conditions i?+ = R-Jr on S/j, with Jr described below. 

(3) R{z) =1 + 0{l/z) as z-^oo. 

The jump matrix J/j for \z\ = pn~^l'^ is by ()6.57p and ()6.56p 

Jr{z)=[i + Z{-z) + ^^ 'j~{z)(i + ^^ 

= {l- Z{-z) - ^) + Z{-^) + Oinr^l^^)) + ^) = / + 0(n-Vi2)^ 

where we also use the fact that Z{z) and Z^^jz are bounded for \z\ = pn~^/^. 

The transformation (|6.57p does not change the jump matrices on the other parts of in 
an essential way. Hence Jr tends to the identity matrix on these parts as well, with a rate of 
convergence that is the same as that for Jy. 

We have now achieved the goal of the steepest descent analysis. R[z) tends to the identity 
matrix as z — t- oo and the jump matrices for R tend to the identity matrix as n — t- c«, both 
uniformly on S/j and in L^(S/j). By standard arguments, see [15] and in particular [9j for the 
case of a moving contour, this leads to the conclusion of our steepest descent analysis 

as n — )• oo, uniformly for z € C \ Tir. 



6.8 Proof of Theorem [M] 

The idea of the proof is similar to that of the proof of Theorem 12.61 i.e. we write the expression 
for the correlation kernel in terms of R instead of Y by unfolding all transformations performed 
in the steepest descent analysis. Since the first few transformations are the same as in the 
regular cases, we see from the arguments in Section [5.81 that if < a;,y < c, 

2m[x — y) 

X (e-"^i.+(^) e-"^2,+W q o)'^. (6.59) 

Moreover, for x,y G D (O, pn^^/'^) , we have by ()6.54p and ()5.60p 

S+{z) = R+{z)S^^\z) = fi+(z)Pl°V2)^diag(e''"*/2^e-^"^/2^e'^"*/2^e-''"*/2^, z = x,y, 
where J is given in (j5.6ip . It then follows that 

„.u/2 -u/2 

2-Ki[x — y) 

X (^p[°\-y)y^ R-\y)R{x)P^°\-x) {O e-"^2.+ W-^-Wa g-nAi,+ W+z^^a)^ . 
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By (inail, IKm . and (fO^D this becomes 

e"(^-^)K„(x,y) = ^- -(0 _^-uni/2\ 

2m{x — y) 

X Ar_ (^n2/3(_y)i/2. n2/3s(_y), ni/3t(_y)) E-^{-y)R-\y)R{x)En,- {-x) 

X N^(n^'^{-xfl'^;ri^-{-x),r2,-{-x),r?l^s{-x),n^l'h{-x)^{<d e^'^^^/^ e'^^/'^ f . 

(6.60) 

Now we scale x and y with n such that 

and y = (6.61) 



^4/3 n^/"' 

where n, t> > 0. Then for large n, x and y belong to the disk D (O, pn^^/^) , so that (|6.6U|) holds. 
We want to take the limit as n — )• oo. Note that under these conditions 



lim e"^^'-^) = 1, 



and by 



J = 1,2, 



n2/3s(z) ^ l(a2 - 56), 
n^/^t{z) 2a, 

as n — 7- oo and z = x,y. Furthermore, it follows from (j6.58p and Cauchy's formula that 

R~\y){R{y) - R{x)) = O = O {n-^y^') , (6.62) 

as n — )• oo where the constant is uniform for n, v in compact subsets of M. Then also 



'R-Hy)R(.x) =(l + Z{-y) + ^) R-\y)R{x) + Z{-x) + ^) ' 

/ + Z{-y) + ^ j ii-i(y)i?(x) {^I - Z{-x) - = I + 0{n-^), (6.63) 

as n — )• oo, where the constant is again uniform for u,v in compact subsets of M. Here we used 
(j6.56p to invert the rightmost matrix. To prove the last equality in (j6.63p . note that the matrix 
function Z{z) — which is analytic by Lemma l6.12f a). can be written as a power series in the 
variable n^/^z with coefficients having a limit for n — )• oo, thanks to (|6.53p and ()6.48p . Applying 
this with z = x,y given in (|6.6ip we get the claimed 0(n~^) estimate in (j6.63p . In fact, the 
same reasoning yields the following more precise version of ()6.63p . 

R-^{y)R{x) = I + P+Oin'^) + 0{n'^)Q+ + o{n-^), n ^ oo, (6.64) 

where the matrices P+, (5+ originate from (j6.53p . (We could also write P_, Q- instead.) 
Next we estimate the factor En{z) given in (j6.50p . We claim that the transformed matrix 

^„(z):=i?„(z)diag((^ ^))diag(C-^/^C'/^C-^/^C^/4), C = n'/'z, (6.65) 

is analytic near z = 0. Indeed it has no jumps, by virtue of Lemma [621 moreover it behaves 
as 0{z~^^^) as z — >• so there is no pole at z = 0. 
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Lemma 6.14. We have 

lim E-\-y)R-\y)Rix)Eni-x) = I, (6.66) 

n— >oo 

uniformly for u,v in compact subsets o/M. 

Proof. In the proof below all the O and o terms will be uniform for u, v in compact subsets of 
R. We start by writing 

E-\-y)R-\y)R{x)En{-x) = E-\-y)En{-x) 

+ E-\-y){R-\y)R{x) - I)En{-x). (6.67) 

Let us estimate the first term in the right hand side. We have En{z)^^ = O (n^/^), z = x,y, as 
n — )• oo, which is a special case of (j6.69p below. Then by the analyticity of En{z) we obtain 

E-\y)(En{y) - En[x)) = 0{{x- y)n) = O (n-^s) , (6.68) 

as n — )■ oo. Consequently the first term in the right hand side of (j6.67p goes to the identity 
matrix for n — ?■ oo. 

Next we estimate the second term in the right hand side of (|6.67p . On account of (|6.6ip . 
and ^^TM we have 

En{z) = P+0 + , En{z)-^ = O (n^/^) Q+ + o {n^'^'^ , z = x,y, (6.69) 

as n — )• oo. From (j6.69p and (|6.64p we see that 

E-\-y){R-\y)R{x) - I)Er,i-x) 

= O {Q+P+0{n-^)P+ + Q+0{n-^)Q+P+^ O (n^l^') + o(l) = o(l), 

for n — )• oo, where the second equality follows from the orthogonality relation = in 

()6.25p . Hence the second term in the right hand side of ()6.67p goes to zero for n — >• oo. □ 

To use the above lemma, we should first express the matrix En in (|6.60p in terms of its 
transformed counterpart £"„ in (j6.65p . This substitution releases an extra factor which multiplies 
from the left the matrix N in ()6.60p . By combining this with the above estimates we find 

lim 4ji^J^,^)=f-^(0 e-/2 _e— /2) 

xiv(^-i^;^/2.2^2, i(a2-56),2a)"^7V(^-iu^/2;2,2, i(a2-56),2a) (0 e'^^^/^ g^^^a^^^ 

(6.70) 

with 

N := diag(cv^c-^/^c^/^c'^/^)diag , -;)) N (-^c^/^) , 

for C = Equivalently, by ([02]) . (lOTD and (|230D . 

N = diag(l, -i, i, l)M(C) diag (l, i, 
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Inserting this in ()6.70p we get 



lim J-i^J^,^)= ^ ^ (0 e— V2 

xM(t;;i(a2-56),2a,i/)"^]w(u;i(a^-56),2a,i/) (0 -e''^*/^ _^g-;^W2)^^ 

where we recah that M(z) = M (z; s, t, u) is defined in ()2.30p with M(C) = M(C; s, t, v) denoting 
the unique solution to RH problem 12.71 with parameters given in (j2.32p . The above kernel is 
clearly equal to the right hand side of ()2.33p . This completes the proof of Theorem 12.81 
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